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ABSTRACT 


This report presents the results of an eight-month study of signal 
processing techniques applicable to the Geodetic Satellite Altimeter program. 

The first subject treated is the analysis of random errors in the altitude 
measurement process which arise from signal fluctuations and receiver noise. 
Results are presented based on both theoretical analyses and computer simulation 
of the altimeter concept. Characteristics of the electromagnetic energy scattered 
from the ocean-surface are then discussed from the standpoint of identifying 
statistical properties of the altimeter signal and for identifying measurement 
biases that may arise in the scattering process. The report concludes with 
a discussion of the presently known oceanographic factors pertaining to the 
scattering problem. 
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SECTION 1 


INTRODUCTION AND SUMMARY OF RESULTS 


L.i INTRODUCTION 

This report presents the results of a study of radar signal processing 
methods for a geodetic satellite altimeter. The two altimeter precision 
requirements considered are: 1) maximum radar system errors of one to two 

meters as required for the GEOS-C program and 2) errors limited to a fraction 
of a meter as required for Sea-Sat A program. Although the emphasis of this work 
is on signal processing, it has been necessary to devote a considerable portion 
of the study effort to. oceanographic and electromagnetic scattering considerations. 

Section 2 deals with the analysis of errors which arise from measurement 
noise and signal fluctuations in the radar implementation. Since these error 
sources are unavoidable in the system, signal processing conditions are discussed 
which reduce these errors to acceptable values. Results are presented based on 
a theoretical analysis and on computer simulations of the altimeter system. 

Section 3 presents a detailed description of the mathematical techniques 
used to simulate radar scattering from the sea surface and to model the radar 
altimeter functions. 

Section 4 summarizes the electromagnetic scattering work performed during 
the study. This subject is of central concern for two reasons: 1) the analysis 
of radar system errors requires accurate characterization of the scattered 
signal, and 2) the identification and compensation of any measurement bias 
arising in the scattering process requires a thorough understanding of the 
underlying physical mechanisms. The two outstanding problems in these categories 
are the modeling of wave height effects in the transient region of the altimeter 
signal and the sensitivity of backscattered power to ocean surface conditions. 

For an idealized ocean surface, e.g. isotropic, Gaussian height distributions, 
the problem has been solved. Section 3 of this report considers the effect of 
more realistic assumptions. 

The mathematical results in the previous section require a number of 
assumptions regarding the ocean surface features. The work reported in 
Section 5 represents a survey of oceanographic literature pertinent to the 
electromagnetic scattering problem. The principal topics considered in this 



Section are the relationship between mean-square slope of the ocean wave 
structure and surface wind, and the spectral features of the ocean. 

1.2 CONCLUSIONS AMD RECOMMENDATIONS 

The system analyses conducted during the course of this study indicate 
that the random errors due to signal fluctuation and receiver noise will result 
in an altimeter precision on the order of one-half meter for a system within the 
following characteristics: 

Single pulse signal-to-noise ratio 10 to 20 db 

Pulse length 50 to 100 ns 

Sampling rate 1000 per second 

These characteristics are well within the state-of-the-art in radar design. 

An altimeter can be designed to meet accuracy objectives of the GEOS-C program 
but particular attention must be given to long-term drift problems. 

There are a number of unknowns in the design of an altimeter with an 
accuracy of a fraction of a meter. One of the most important questions aside 
from "sea state" bias is the pulse-to— pulse correlation of the altimeter signal. 
This limits the attainable altimeter accuracy per unit time. Radar reflections 
from the sea have never been accurately measured under satellite conditions and 
the GEOS-C satellite is probably the best approach to obtaining this information. 

Although the GEOS-C performance can be realized with a conventional pulsed, 
split-gate, or threshold signal processor, the Sea-Sat A equipment will require 
more sophistication. The theoretically computed random error for a 50 ns pulse 
length and for 1000 samples per second is 20 cm. This error can in theory be 
further reduced through use of shorter pulses, faster pulse rise-time, or more 
elaborate transmitter waveforms, even if the pessimistic assumption of a one- 
millisecond signal correlation time is found to apply to satellite data. The 
results given in this report indicate that systematic errors and equipmental 
biases will constitute the largest instrumentation error in the Sea-Sat A 
concept. These non-random errors can arise from effects such as: 1) mean-value 

shifts in the altitude data as a function of signal statistics or signal-to-noise 
ratio, 2) environmental or temporal drift characteristics of -the satellite 
equipment, or 3) unrecognized processor non-linearities. Because of the severity 
of these problems it is recommended that a wide range of radar techniques be 
investigated for Sea-Sat A. It is further recommended that future radar altimeter 
research emphasize the Sea-Sat A requirements, since added knowledge of problem 
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areas and techniques required for the Sea-Sat A system may lead to a more 
evolutionary concept for GEOS-C. Until the first satellite altimeter is in 
operation, many elements of the altimeter function will remain speculative. 

In regard to electromagnetic scattering, it is found that the radar 
cross-section O q can be related to mean squared slopes of an isotropic sea 
surface. For non-isotropic ocean surface conditions the relationship is more 
complex. Derivation of the functional relationships through a theoretical 
electromagnetic approach appears unrealizable at this time because of the 
extreme difficulty in obtaining accurate high frequency ocean-wave data. 
Equivalently, the ocean surface autocorrelation function cannot be measured 
with the required accuracy using existing methods. An empirical approach is 
therefore recommended for obtaining normal incidence data in which actual radar 
data is correlated with ground truth information under varying sea surface and 
meteorological conditions. Because of the normal incidence geometry problems 
and altitude limitation associated with conventional aircraft measurements, 
extraction of such data from the GEOS-C experiment is strongly recommended. 

For the investigation of "sea state" effects on the altimeter signal, 
acquisition of near-surface (short pulse) radar and laser profilometer data is 
recommended. Such data would constitute a basis from which to assess the effects 
of the approximations and assumptions in the electromagnetic models of sea-state 
bias . 

The principal conclusion of the oceanographic study is that mathematical 
arguments require the two-dimensional power spectrum to exhibit 180° symmetry. 

At present, there is no single technique which will provide a two-dimensional 
spectrum of the accuracy and spatial resolution needed for electromagnetic 
scattering investigations. 
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SECTION 2 


RADAR SYSTEM STUDY 


2.1 BACKGROUND 

This section presents the results of the study pertaining to radar signal 
processing. A number of important error sources have been investigated during 
the course of the program. These include: errors arising from the signal 

fluctuations inherent in planetary or ocean scattering*, errors resulting from 
the limited number of samples available per unit time, and errors caused by 
thermal noise and processor non-linearities. The importance of these errors is 
assessed relative to the two altimeter precision categories and to techniques 
for minimizing and/or compensating for these errors. 

This section is organized as follows: As a means of establishing concepts 

and nomenclature, a general description of the altimeter techniques under 
consideration is given. Computed waveforms are shown to clarify the concepts. 
This discussion is followed by a presentation of the principal results of the 
radar altimeter study. Error sources and parametric effects are investigated 
using theoretical results and computer simulations. The section concludes with 
a consideration of general system characteristics and a review of potential 
alterations to the radar system. Computational aspects of the simulation and 
radar characteristics which are heavily influenced by either electromagnetic 
or oceanographic considerations are considered in later sections. 

2.2 DESCRIPTION OF THE -PRESENT ALTIMETER CONCEPT 

A number of organizations have considered the problems of measurement of 
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satellite altitudes to the precision required in the geodetic investigation . 
The more conventional system characteristics such as transmitted waveform, power 
level, sensitivity, bandwidth, and antenna gain have been covered in the cited 
references and will not be discussed here. As presently envisioned, the first 
generation altimeter will consist of an X-band pulsed radar with provisions 
for measurement of time-of -arrival of the received signal. The development of 


* The term "self noise" is commonly used in radar astronomy to describe 
this effect. ’ 
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such a system differs in several important areas from the design of conventional 
radar systems. Attainment of the desired accuracies will require optimum signal 
processing and timing techniques, and knowledge of the effects of oceanographic 
features on the scattered signal. 

Some of the technical problems that must be considered in developing 
satellite equipment are: 1) the limited space, power, and weight available, 

2) the nature of ocean backscattered signal, and 3) satellite dynamics. As 
discussed more fully in Section 4, return from the ocean surface arises from 
many discrete regions. The signal is dispersed in both time and frequency. For 
the assumed satellite conditions, a 50 nanosecond (ns) transmitted pulse will 
give rise to a received signal with a time spread of 2. 6 microseconds for a 
3 degree antenna beamwidth. The signal will fluctuate in amplitude with 
characteristics similar to signals received over a rapidly fading channel. 

These signal characteristics can be seen by an examination of the receiver wave- 
forms, shown in Fig. 2.1, which were obtained as a by-product of the simulation 
study. Figure 2.1 shows ten typical received waveforms corresponding to the trans- 
mission of a 50 ns rectangular pulse scattered from the ocean's surface and 

received by a very wideband receiver. Waveforms with the vertical scale labeled 

2 

E represent the output of a linear envelope detector and those labeled E 
represent the output of a square-law detector. The horizontal scale shows 
relative time in nanoseconds, with 25 ns corresponding to the instant at which 
one-half the pulse envelope arrives on the sea surface. These computed waveforms 
correspond to independent samples. For a simulated radar inter-pulse period 
less than the correlation time, the waveforms would show evolutionary changes. 
Referring to Fig. 2.1, except for the transient region in which the entire 
pulse is not incident on the surface, the sea return signals are much like 
samples of receiver noise. 

Simulated waveforms are shown in Fig. 2.2 corresponding to a Gaussian 
shaped pulse that is 50 ns wide at the e ^ points. The horizontal scale shown 
in this case is based upon the center of the pulse arriving at the ocean surface 
at 51 ns. The same distribution of reflecting area was assumed in the computa- 
tions used in Figs. 2.1 and 2.2; the Gaussian returns therefore appear as 
smoothed versions of the rectangular data. The Gaussian envelope is considered 
indicative of matched filter IF characteristics. Note that the Gaussian results 
demonstrate that filtered waveforms need not be sampled rapidly for good 
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Fig. 2.1. Typical Simulated Results of Square-Law and Linear 
Detector Waveforms for a 50 ns Rectangular Pulse, SNR = «. 












reconstruction of the waveform. Also, it is noted that for a system containing an 
AGC response- time on the order of one second, some of the received waveforms 
will contain very little average energy over a 200 ns time interval. Case 5 of 
Fig. 2.2 demonstrates the fact that the instantaneous altitude error can be 
several times the pulse length for a matched filter IF regardless of the measure- 
ment scheme used. Some of the waveforms shown in Fig. 2.1 exhibit a saturation 
effect. This is not a result of the simulation program; is is simply due to the 
scale size used in the figure. 

Figure 2.3 shows the effect of receiver noise on the altimeter waveforms. 

The same signal characteristics were used as in Fig. 2.2, and a comparison of 
Figs. 2.2 and 2.3 shows the effects of a 10 db signal-to-noise ratio (SNR) on 
the waveforms. The fine structure in the noise is due to the finite slope 
of the noise spectrum with frequency. 

The double-delay differencing type of altitude processor consists of two-stage 
signal differencing with a delayed and inverted replica of the original signal^. 

The sample waveforms involved in the double-delay differencing operation are shown 
in Fig. 2.4 (the signal waveform previously shown as Case 1, Fig. 2.3 was used). 
Figure 2.5 shows typical waveforms at the output of the double-delay differencer 
for ten individual cases (for the noise free Gaussian signal) . A comparison of 
Figs. 2.4 and 2.5 shows that individual cases depart drastically from the results 
that would be obtained by using an idealized ramp signal. The most significant 
feature in Figs. 2.4 and 2.5 is that multiple zero-crossings are present. 

Because of these ambiguities, the double-delay technique would not be suitable 
in the altimeter without the addition of a threshold circuit or a closed-loop 
impl ementation . 

Additional examples of simulated waveforms are contained in Appendix A. 

The quantitative results of the radar system study are discussed in the next 
section. 

2.3 DISCUSSION OF ERRORS IN THE ALTITUDE MEASUREMENT PROCESS 

The two principal methods of altitude extraction examined in this report are: 

1) double-delay differencing and 2) thresholding. The investigation was limited 

to these techniques, because the General Electric Company had concurrently 
• * 4 

investigated the split-gate technique . The results given below apply to both 
open and closed-loop systems. 
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Fig. 2.3. Typical Simulated Results of Square Law and Linear Detector Waveforms for a 50ns Gaussian Pulse 

SNR= 10 db 











Histograms of the individual altitude errors at the output of a double- 
delay differencer were obtained through digital simulation for fifty independent 
cases. A 50 ns Gaussian shaped pulse was used; the delay interval for the 
double-delay difference was 50 ns. Figures 2.6 and 2.7 show the histogram data 
for signal-to-noise levels of <», 30, 20, and 10 db. The standard deviation o 
and mean crossings y derived from the histogram data are shown in Table I. 

The altitude error to be expected for a 1000 sample average also was 

obtained from the histogram standard deviation o through the equation 

s 

a = 1.62 — a 
s 2v'n s 

where n is the number of samples, c is the propagation velocity, t is the unit 

of time of a , and the factor of 1.62 is the conversion of pulse length used in 
S -1 

the simulation (e points) to half-power pulse length. 


Table 2-1 

Standard Deviation and Mean of Histogram Data 


SNR 

Std. Dev. 
(Og) in ns 
50 cases 

o s in 
meters 
1000 cases 

Mean 

(y) in ns 

CO 

18 

.14 

94.6 

30 db 

18 

.14 

96 

20 db 

17.2 

.13 

89 

10 db 

18.4 

.14 

68 


Table I shows the altitude variance to be largely independent of receiver 

noise. The indicated altitude uncertainty a converts to a (two-way) altitude 

s 

error of approximately 0.14 meters for 1000 samples. 

A theoretical analysis of the altitude error characteristic is given in 
Appendix B. This analysis assumes' that correlation exists within the detected 
altimeter signal (i.e. matched filter IF) and that a double-delay differencer 
type of altitude measurement is used. The altitude error 0 in meters is found 

Si 

to be 
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Fig. 2.6. Histograms of Pulse-by-Fulse Altitude 
Jitter, 50 Computed Cases 
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Fig. 2,7. Histograms of Pulse-by Pulse Altitude 
Jitter, 50 Computed Cases 




o = -J^\^68 + 9.5 (SNE)" 1 + 7,9 (SNR)“ 2 

a t- Y 

/n 

where T is the pulse 'length in .nanoseconds (at-. the half -power point) and n is 
the number of independent samples. A comparison of theoretically computed and 
simulated altimeter errors (arising from signal fluctuations and receiver noise) 
is shown in Fig. 2.8. The two results are in good agreement. ‘ 

The mean-value figures given in Table I show considerable change with 
noise level. This fluctuation may be caused by two effects: 1) variation due 

to the limited sample size present and 2) variations caused by any differences 
in the spectral characteristics of the signal and receiver noise. To determine 
which effect was dominant, a series of computer runs were made under varying 
conditions. Figure 2.9A shows the histogram results obtained for a 20 db SNR 
and for new sets of random numbers used in the signal and noise simulations. 

Figure 2.9B shows histogram results for 20 db SNR using a narrower noise bandwidth. 
These results are summarized in Table II. 


Table 2-2 

Effect of Simulation Conditions on the Values Shown in Table I 

(SNR = 20 db) 


Computation 

Std. Dev. 
(a s ) in ns 
50 cases 

S a in 
meters 
1000 cases 

Mean 

(p) in ns 

Original 

1-7.2 

.13 

89 

New random numbers 

22.8 

.18 

96 

Reduced noise bandwidth 

17.9 

.14 

91 


The mean and variance are both seen to be more dependent on initial conditions 
than on noise bandwidth. A much larger sample size (~1000) is therefore required 
to provide acceptable accuracy in the histogram data. Computations of such 
magnitude were outside the scope of this contract and can be performed more 
efficiently using analog techniques. 

The fact that the mean crossing point is sensitive to receiver noise 
characteristics may be seen from exaniming the theoretical relationship for the 
probability of a zero-crossing p(x) . For Gaussian processes which start at t = 0, 
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Fig. 2.8. Comparison of Theoretically Computed and 
Simulated Altitude Errors Versus SNR 
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Fig. 2.9. Histogram of Pulse-By-Pulse Altitude 
Jitter, 50 Computed Cases 
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p(x) is given by , 

p(t) .iJmSsLsJMl 

' ’ E(o) 

in which R(x) is the correlation coefficient. For additive (uncorrelated) 
signal and noise this expression becomes 

. R (t) + R (t) " 
s n 

' R (o) + R (o) 
s n 


P( T ) = ~ 

7f 


v 


where the subscripts s and n refer to signal and noise. Assuming exponential 
correlation functions of the form 

-at -at 

s . n 

e and e 

-a t -a t- 
S it 

for a^ > a^, the numerator of the square root term e + e decays faster 

than for a^ = and the zero crossing probability occurs earlier in time. 

Any shift in the mean value of the altitude indication constitutes an . 

important source of error in the altimeter. The above results, although 

qualitative because of the sample-size limitations, indicate that a bias will 

exist that is a function of SNR. A similar effect has been noted in the split- 

4 

gate tracker when a dc restore circuit is not used . 

Neglecting wave-height effects, the integration of a large number of radar 
returns from the sea-surface should produce a "ramp-like" signal since the area 
exposed during the early portions of the pulse illumination will increase linearly. 
A presentation that is proportional to signal power (i.e. a square-law detector) 
will therefore exhibit this linear dependency. These computed waveforms for both 
linear and square-law detectors are shown in Figs. 2.10-2.12. Figure 2.10 shows 
the waveform obtained by averaging the 50 rectangular pulse cases (see Fig. 2.1). 
Confidence bounds may be estimated for these data by noting that for a Rayleigh 
distribution (linear detector case) the variance is equal to Cl - times the 
mean value p, and for an exponential distribution (square-law detector) the 
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Fig. 2.11. Averaged Waveform for Gaussian Pulse (Original Computation)^ 50 Cases 
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variance is twice the square of the mean value. Therefore, standard deviation 
for the 50 sample waveform is 


a 




Rayleigh 



Exponential. 


These one-sigma hounds are shown as dotted lines in Fig. 2.10. Figures 2.11 
and 2.12 show the average waveforms for two simulation cases of 50 samples each, 
using Gaussian shaped pulses. 

Estimates of the errors to he encounterd in altimeter operation have been 
made using mean value waveforms as in Fig. 2.12. One may argue that the averaging 
process can be performed prior to double-delay differencing and certainly the 
mathematical and conceptual simplifications are appealing. The result of double- 
delay differencing of an average of 50 waveforms is shown in Fig. 2.13. The 
zero crossing for this average waveform occurs at 99 nanoseconds. The results 
of averaging single pulse errors (Table I) indicate that the mean zero crossing 
occurred at 94.6 nanoseconds. Significant differences can therefore exist 
between the two computations, and analyses based upon an equivalence between the 
two should be treated with caution. 

2.4 EVALUATION OF THRESHOLD TECHNIQUES 

Simulation results are shown in Figs. 2.14 and 2.15 for a threshold type of 
altitude measurement. The threshold levels used in Figs. 2.14 and 2.15 were, 
respectively, 50 percent and 33 percent of the maximum value of the average 
signal. These results contain several interesting features. The 33 percent 
threshold level (Fig. 2.15) results in an altitude standard deviation of 18 ns 
(50 cases) , which is exactly the result obtained using a double-delay differencer 
(c.f. Fig. 2.6A). The 50 percent threshold was found to produce a larger standard 
deviation (25 ns) because of occasional values with large time deviation (the 
tail shown in the histogram) . The fact that some of the threshold crossings occur 
long after the leading edge can be seen by examination of the waveform data in 
Case 5 of Fig. 2.2. This particular waveform crosses about 140 ns late. The data 
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Fig. 2.14, Histogram of Pulse-By-Pulse Altitude 
Jitter for 50% Threshold 
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Fig. 2,15. Histogram of Pulse-By-Pulse Altitude 
Jitter for 33% Threshold 
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suggest that the 50 percent threshold performance could be improved by either 
receiver AGC action or by discarding or limiting the very late crossings. This 
last technique is somewhat akin to the limiting behavior of the time-discriminator 
used in split-gate systems. Use of AGC would require a sufficiently high prf so 
that AGC action could be derived on one pulse and applied to the next (correlated) 
pulse. 

The simulation results show that both the 33 percent and 50 percent thresholds 
result in mean-values that are very close to the ji priori values (levels selected 
from average waveforms such as 3?ig. 2.12). In view of the dependence upon SNR 
of the mean value encountered with the double-delay circuit, the threshold data 
suggest that a threshold technique may be less susceptible to this type of bias. 

In the future, threshold simulations will be conducted to determine the effect 
of SNR upon the mean value of the output. 

Before leaving the threshold technique, mention should be made of the false- 
alarm problem and the open-loop altitude measurements. In order to reduce the 
false alarm rate, the threshold technique would require a much higher single pulse 
SNR than would the split-gate technique. Possible solutions to this problem, 
within the power constraint of a TWT transmitter, are: 1) use of pulse compression 

or 2) use of surface wave transversal filters . Taking into account pulse 
stretching due to the antenna pattern, the transversal filter method could be 
implemented by transmitting several pulses spaced three microseconds apart. The 
threshold circuit is inherently suited to open-loop filtering, whereas the split- 
gate and double-delay circuits are basically closed-loop sensors. 

2.5 EFFECT OF VIDEO NON-LINEARITY 

In order to test the effect of non-linearity on a double-delay differencing 
form of altitude measurement, a video transfer function of the type 

E = E . - E? 

out xn 25 xn 

was simulated. Such non-linearity was assumed to be unknown (e.g,, due to drift) 
in the altimeter. This function represents a quadratic non-linearity in the 
processor. Results of this simulation are given in Fig. 2.16. The variance is 
seen to be unchanged and the point of mean zero crossing shifted 1.4 ns compared 
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to the data given in Fig. 2.6A. The assumed non-linearity represents a one 
percent departure at the 50 percent point on the mean waveform. The resultant 
shift is therefore somewhat larger than mean-value estimates would predict. 

For altimeter accuracies of 1 - 2 meters, non-linearities of this magnitude are 
not highly significant. For the Sea-Sat A, such non-linearities will be important. 

2.6 PRF AMD MISCELLANEOUS CONSIDERATIONS 

Correlation properties of the altimeter signal may be estimated through 
computation of the Doppler spread of the scattered signal. Previously published 
estimates indicate that the Doppler spectrum will be in the neighborhood of 
1600 Hz\ These figures are based on the signal characteristics at the break 
between the ramp and the plateau of the mean waveform, for a 50 ns pulse. 

For the proposed GEOS-C orbit, a maximum vertical velocity component of 
about 700 ft. /sec. will be present and will result in a (maximum) Doppler shift 
of approximately 10 KHz. Estimates of the signal correlation properties derived 
from the Van Gittert-Zernike Theorem"^,, place a value of up to 10 KHz on the 
signal fluctuation. For the purpose of this study, a conservative value of 1000 Hz 
was used. 

Measurement of the vertical Doppler component would facilitate altitude 

extraction with elliptical orbits, .since rate prediction could be used in the 

signal processor. Alternatively, it might be possible to difference the smoothed 

(e.g. one second) estimates of altitude to obtain prediction information. The 

third alternative would use a ground command and control function for this purpose. 

These techniques could minimize the problems associated with bandwidth and aquisi- 

4 

tion in a closed-loop tracker . The ground control technique appears preferable 
from the standpoint of minimizing satellite equipment complexity. 

Appendix C presents the results of a brief survey of prior work in optimal 
processing. The optimal video processor requires an impulse response that is 
related to the time-inverse of the sea-surface scattering function. From an 
engineering standpoint, it should be preferable to utilize only that part of 
the scattering function in the neighborhood of the ramp region. Accordingly, the 
optimum processor closely resembles an integrate-and-dump circuit, such as is in- 
volved in the split-gate tracker. It is therefore doubtful that any significant 

Suggested by G. Bush of Johns Hopkins University, Applied Physics Lab. 
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advantages would accrue from synthesizing the .theoretically optimum processor 
encountered in the survey. Adaptive filtering.should.be investigated. 

2.7 SUMMARY 

A comparison of the altitude uncertainty for three sensing techniques is 
shown in Table III. 


Table 2-3 

Comparison of Altitude Errors 
50 ns Pulse 1000 Samples 


Sensor 

Alt. Error 
(meters) 

Double-delay (20 db SNR) + 

.155 

Split-gate (20 db SNR)* 

.22 

Threshold 50% (<» SNR) 

.25 

Threshold 33% (20 db SNR) 

.14 

Theoretical (20 db SNR) 

.20 


Average of first two entries shown in Table II 

* Scaled from data for a 100 ns pulse system with dc restorer, given in 
Ref. 4. 


The results summarized in Table III indicate that any of the three sensing 
techniques considered are usable for the GEOS-C altimeter. (The data for the 
split-gate technique, shown in Table III, was not derived from the simulation 
and the figures shown are not directly comparable.) Because of the more extensive 
hardware experience and greater mean-value stability of the split-gate technique, 
it is recommended for closed-loop implementations. The threshold technique merits 
further work, especially if studies 6f Sea-Sat A requirements show non-analog 
techniques to be useful. From the processor standpoint, the main problems with 
Sea-Sat A are due to the mean-value shifts (e.g., non-linearity,, systematic 
changes with SNR, environmental effects, and temporal variations) . More informa- 
tion is needed regarding altimeter signal correlation properties (pulse-by-pulse) 
to determine the best solution to the signal fluctuation error problem with Sea-Sat A. 

The theoretical analysis in this report shows that self-noise is a larger 
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2 rror source than receiver noise, for a SNR greater than 10 db (see Fig. 2.8). 
Ehis analysis correponds to the use of matched (non-adaptive) IF bandwidths and 
should be extended to treat pulse rise-time as a variable. A wider bandwidth 
per se will reduce self -noise errors. Altitude variances for a threshold 

system using rectangular pulse data (given in Appendix B) have been found to 
tesult in one-half the variance obtained for Gaussian shaped pulses. 

In summary 5 the altimeter accuracy requirements for GEOS-C are within the 
state-of-the-art in radar performance. It is recommended that future studies 
amphasize Sea-Sat A problems and the development of sampling techniques for use 
with GEOS-C. The availability of information such as coherence time from the 
GEOS-C experiment would be of great value to the Sea-Sat A program. 
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SECTION' 3 


ALTIMETER SIMULATION MODEL 


3.i GENERAL DISCUSSION 

The preceding section of this report presented results of theoretical 
analyses and computer simulation studies of the altimeter radar concept. 

The Monte Carlo technique used to generate the simulation results is described 
in detail in this section. The simulation approach was indispensable in several 
areas of the altimeter study because: 1) the mathematical complexity of the 
phenomena under investigation did not permit closed-form solution, or 2) the 
assumptions and approximations contained in available theoretical results led to 
limited confidence in the solutions so obtained. The statistical level-crossing 
problem and the analysis of time-varying signal statistics are examples of radar 
system functions that have not been rigorously solved. For non-Gaussian 
statistics there is a paucity of even approximate relationships in the literature. 

The overall simulation program is shown in block-diagram form in Fig. 3.1. 

It consists of the computation of the ocean-scattered pulsed signal characteristics 
and the modifications resulting from the receiver bandpass characteristics, 
additive thermal noise, detector and other non-linear characteristics, and the 
altitude measurement process. Waveforms and statistical compilations are then 
obtained from the output data. The salient operations shown in Fig. 3.1 are 
described in detail below. 

The electromagnetic scattering process for the ocean's surface may be 
visualised as due to specular regions superimposed on the gross wave structure. 
These specular regions consist of areas with proper orientations and curvature 
relative to^ the wavelength of the illuminating radiation. The scattering process 
is governed by the high frequency region of the ocean-wave spectrum. As discussed 
in Section 4, backs cattered energy at normal incidence is a function of the mean- 
square slope of the ocean’s surface. Since the backscattered signal from the 
ocean's surface is due to contributions from many individual scattering elements 
within the radar beamwidth, the well-known Rayleigh distribution of amplitudes 
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Fig. 3.1. Block-Diagram of Computer Simulation 





results''". The principal reasons for the wide acceptance of this distribution 
are: 1) the fact that observed .returns provide a reasonably good fit to Rayleigh 
statistics, and 2) the plausibility of the physical model (i.e. the specular 
model) . 

The simulation of ocean scattering is based upon the characteristics for 

backs cattered power at normal incidence. For plane wave illumination at normal 

-2 

incidence, the average backscattered power E g at a distance R q from the scattering 
surface is 2 



s 


M 2 e 2 

4irR 2 



(AREA). 


( 1 ) 


where 

K is the Fresnel reflection coefficient 
is the incident electric field intensity 

2 

h Q is the variance of surface displacement about mean-sea-level, and 
Y is related to the sea-surface correlation length. 

This equation indicates that for a given radar wavelength and a time-invariant 
sea, the backscattered power is proportional to the illuminated area. Since 
the area illuminated by a finite antenna beamwidth impinging upon a spherical 
earth is a function of time, the backscattered power is also proportional to time. 

The area illuminated by a satellite-borne radar pointed toward the nadir is 
■u 2 

given by 


A(t) - ir a c t for t < t 

“ P 


where 


a is the satellite altitude 
c is the velocity of light, and 
t is the transmitted pulse length. 


For t > t , 
P’ 


the radius of the leading edge is given by 


r 


L 


(a c t) 


1/2 


( 2 ) 
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and the radius of the trailing edge is given by 


r t * la c(t - t p )] 1/2 . 


The projected illuminated area then is 


A(t) 



7T a c t for t < t 

- p 


(.3) 


(4) 


and 


A(t) = 




ir a c t for t > t . 
P P 


(5) 


Thus, the projected area increases linearly for the length of a transmitted 
pulse, after which it remains relatively constant until other effects dominate. 

3.2 MATHEMATICAL DESCRIPTION 

In order to utilize the above model of radar return, macroscopic conditions 
must prevail since the existence of specularly reflecting points is a random 
occurrence. For the geodetic satellite geometry, examination of the scattered 
signal in one nanosecond increments corresponds to area increments of one square 
kilometer. An area of this size should contain at least five independent scatters 
and the Rayleigh distribution will match the expected signal statistics. 

The Rayleigh amplitudes are simulated for each one-kilometer area by generating 
numbers and y^ from independent normal distributions with zero mean and unit 
variance. The Rayleigh signals e(n)^ at the output of a linear detector for 
rectangular pulse illumination are 


* w j ‘ 




1/2 


n < t 


( 6 ) 


where the n subscript indicates the time index in one-nanosecond steps for the 
jth sample function. For a rectangular pulse, of length n = 50, the partial sums 
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(n £ '50) of e(n) . represent the transient region of the scattered pulse. For 
example, the partial sum for. the first interval e(l).. corresponds to the aggregate 
return from the reflecting areas exposed during the first one nanosecond. 

(The term "time interval" and one nanosecond increment will be used interchange- 
ably in the following discussion, although time can be scaled to an arbitrarily 
larger interval. Smaller time intervals correspond to areas for which the 
ocean surface statistics may not be independent.) For n = 2, e(2)_. corresponds 
to the vectorial magnitude of returns from both area increments. Since the 
scattering areas are assumed uncorrelated, signal power will be additive and 
e(2)_. will increase linearly (in the mean) up to the assumed pulse length (n = 50). 
- - ■ -The simulated waveforms obtained for a rectangular pulse are mainly of 
theoretical interest because of the large receiver bandwidth implied. Figure 2.1 
showed such waveforms for both the linear and square-law detector outputs. In 
order to simulate more realistic conditions, a Gaussian shaped pulse was used in 
the study of radar parameters. This waveform was chosen because the output of 

3 

real networks with near-rectangular input signals resembles a Gaussian waveform . 

A Gaussian envelope f(t) of width W measured at the e points is given by 


f(t) 


and the width of the spectrum A(f) 


A(f) 



■ “ p [- m] 2 

, -I 

at the e points 



for this envelope is 
(i2mft) dt 



(7) 


The detected waveform for a Gaussian shaped pulse can be generated by the 
expression 


e(k)j = 




exp 
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( 8 ) ■ 



in which the time index k begins at k ' -50 and terminates at k = 150. The 
peak of the pulse in this case corresponds to k = 1. 

Equation (8) assumes that the effect of a Gaussian shaped transmitted 
pulse can be combined with the effect of a .receiver containing a Caussian 
shaped impulse response. This assumption is based on the linear model for sea 
return and 'the fact that convolution of two Gaussian functions yields a 
third Gaussian function. The convolution of the form 


and 



(9) 


can be shown to result in the Gaussian function 


x 3 (t) 


exp 


-t 


2 2 
T + T 

2 


( 10 ) 


Assuming post detection or video filtering can be described by an RC type filter 
with corner frequency given by 


f l ~ 2 ttRC 


(ID 


then the inpulse response h(t) for the filter is 


h(t) R(J exp J- RC 


| = 2n-f ^ exp | -2-n-f^t j- 


( 12 ) 


If f^ is equated to a multiple m c*£ the one-sided bandwidth of the detected pulse 
given by equation (7), i.e. 


f 


1 


2m 

' 71 W 
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then h(t) is 


h(t) = ^ exp j- t} . (13) 

Combining the above results, the signal at the output of an ideal 

envelope detector E is 

n 

E(n) i “ V k S 5Q e <»;| “P j- n? <■> - *> J • (14) 

2 

Similarly, the output from a square-law detector E (n) is given by 

s2(n) j = ^ k ^ 0 e j (k) exp I " (n “ k) | • (15) 

For the simulation results given in this report, k varied from -50 to +150 
and W = 50. This corresponds to a 50 ns pulse, the peak of which arrives on -the 
sea-surface at k = 1. Receiver time delays are not included in the formulation. 
Receiver noise is added to the calculations through the equation 



in which c^ and s^ represent in-phase and quadrature noise voltages selected 
from independent distributions. These distributions are correlated in time to 
simulate the effects of receiver predetection bandwidth characteristics. 

Before considering the details of receiver noise generation, the remaining 
aspects of the simulation will be discussed. The satellite altitude information 
is obtained by processing the video signal to extract an estimate of the time 
between signal transmission and receptions. The principal schemes for extracting 
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this estimate are: 1) early gate, late gate, 2) double-delay differencing, and 

3) threshold detection. The double-delay differencing and threshold techniques 
have been implemented thus far in the simulation. 

The threshold function is accomplished by level filtering in the digital 
program. The double-delay differencing operation is implemented through the 
discrete-time equivalent of the impulse response 

h(t) = $(t) - 2<$(t - x) + 6(t - 2x) 

where x is the time lag and 6 is a delta function. In programming this equation, 
the delayed replicas of the signal are combined with the original input through 
computer storage. 

The remaining operations performed in the simulation program consist of 
non-linear transformations and the compilation of statistical descriptions 
(such as histogram data) . The former is used to study the effect of system 
non-linearities on the altimeter data. To date, saturation forms of non-linearities 
of the type 


E* » E - kE 2 
n n n 


(17) 


have been used, in which k is a parameter. 

In order to establish signal-to-noise ratios (SNR) in the simulation, the. 

2 

variance of the signal process a must be known. For a Gaussian shaped signal, 

s 

the autocorrelation function (ACF) of the signal R (x) is 

ss 


E ss ( t) - 




CO 

exp 


{-B 2 (t - x) 2 } exp {-B 2 t 2 } dt 


(18) 


in which k is the spectral density function for a shaping filter realization 
s 

and B = 2/W. This function integrates as 


R Co) 
ss 



d 2 


= a 


s 


C19) 
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and for B = 2/50 .the signal variance is seen to be 


°s ■ k s < 25 > Vf • C20) 

Two techniques have been developed for the digital noise computations. The 
simulation results given in this report are based on the technique to be discussed 
first. The second method, which has not been programmed, was developed to assess 
the sensitivity of the zero crossing probability distributions to the noise process. 

The first technique generates a correlated noise sequence through a recursive 
first-order difference equation which is driven with random noise. The c^. and 
terms in (16) are generated through separate difference equations of the form 

c t “ Vk-i + w k (21) 

in which is an uncorrelated sequence (zero mean) . The quantity A^ is obtained 
by forming the expectation E of c^, i.e. 

E [ C k-l C lJ = E [ C k-l (A k C k-l + w k^ = \ E [ C k-l C k-l3 (22) 

since Etc^-jW^] is zero. Assuming that is a zero mean process, this last result 
can be expressed as an ACE, i.e. 


or, in general 


R k,k-1 W-l.k-l 

a = i 

k " \k " p J k 


(23) 


(24) 


which is the normalized ACF. The sequence c^ has a Gaussian realization with an 
ACF 


p jk 


e - ct|k-j | d s -o|T 0 | 


(25) 
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The variance of the process can be found through the development 

«¥k> - E t c Vk-i + v ( Vk-i + w kd 

“ E C^k-l 5 + E < M k> 

(26) 

= E [4cVl C k-2 + w kfl )2 ] + EC "k> 


If A is time- invariant, this becomes for a unity variance w^. 


E(c k c k ) = 1 + A 2 + A 4 + . . . 


(27) 


which can be summed as 


E< Vk> - 


1 + 


= a. 


1 - A 


N 


A < 1 


(28) 


This equation is used to specify the steady state variance of the c k sequence. 
Notice that the recursion must be initialized from a random number with the 
proper variance if the computations are re-started. 

The correlation factor a in the ACF was obtained through a shaping filter 
analogy. If a shaping filter H(jco) driven by white noise produces an exponential 
ACF, then the output spectrum <j>(w) must be 

r co _ 2 

<Kco) = R(T) e~ ja)T dT = 2 a 2 . (29) 

J -a> co + a 

The transfer function of the shaping filter is readily seen to be 


H(jco) 


1 


.0) , , 

j b 1 

J a 


(31) 
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Note that a is the 3 db bandwidth in radian measure. The equivalent rectangular 
bandwidth is 



(31) 


Thus the H(jw) transfer function has a bandwidth equivalent to an idealized 
rectangular bandpass of (iT/2)f. 

For a comer frequency of 20 MHz, the ACF is 


P ( T ) 


-125 x 10 6 t 

= e 


(32) 


where t is in nanoseconds. 

2 

Analogous to (18) , the signal-to-noise ratio is obtained by defining cr as 


oo 9 

2 _ 1 a _ ka 

°n 2 tt 3 a2 + m2 U 2 

J —CO 


(33) 


Therefore, the signal-to-noise ratio for the simulation is 


2 

a 

SNR = -| = 

a 

n 



(34) 


The second method of noise generation was designed to provide noise correla- 
tion properties that are matched to the signal characteristics. It assumes that 
the IF bandpass characteristics are Gaussian to first-order. If a Gaussian band- 
pass to an Nth order approximation were used, then the normalized ACF (p(x)) 
could be computed for all time lags T. Noise would thus be generated by considering 
not just the last value, but all previously known values in a. joint probability 
distribution of N variables. For a first-order approximation, it was assumed that 
p(t) is known for only a fixed increment t. 

For the bivariate Gaussian distribution 


P(I 1’ I 2 ) 


1 

2ira 2 [l - p 2 (t) ] 


exp 



ij.+ .l 2 - 2l x l 2 p (- t) " 


j 

2« 2 I1 - p 2 (x)] 



(35) 
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where the a variances' are equal, letting = 1^ + AI and integrating over 1^ 
gives 


P(AI,T) 


2ayir(l - p(t)) 


exp 


{■ 


2 


(AI) 

4a 2 [1 - p(t)] 


(36) 


in which P(aI,T) is the probability of a change in I in a time T. This distribu- 

2 

tion is Gaussian with variance 2a [1 - p(t)]. 

The numerical values of the correlation parameter can be established as 
follows. For a Gaussian waveform 


y(t) = exp 



(37) 


in which T is one-half the pulse duration at the e points, the spectral width 
(one sided) is f Q = 1 /(Ttt). Using a transfer function with a squared magnitude 
of 


G(f ) | = exp 


1 -2 


(38) 


the normalized ACF p(x) is found to be 


■ !• 
J — CI 


)(t) = exp 


2f 


exp {-2nifT} df = exp 


| | 


(39) 


Therefore, for 2T = 50 nanoseconds, the correlation parameter is 

2i 

l /„ \ 

p(x) = exp 


-(37) 


(40) 


in nanoseconds. 
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SECTION 4 


REVIEW OF ELECTROMAGNETIC SCATTERING FROM ROUGH SURFACES 
AS RELATED TO THE GEODETIC ALTIMETER 

4. i GENERAL DISCUSSION 

Satellite altimetry, as presently proposed for the GEOS satellite, 
is heavily dependent on characteristics of the backscattered signal 
from the ocean surface. The altimeter techniques under consideration use 
a beamwidth sufficiently broad (2-3 degrees) that with reasonable satellite 
attitude control, strong radar returns will always be reflected from the 
sub-satellite point. The use of pulses in the range of 50-100 nanoseconds 
ensures that at satellite altitudes of 1000 km, the transient portion of 
the return will indeed be backscattered normally from the sea surface. 

At satellite altitudes of 1000 km, the footprint of the initial 

100 nanoseconds of the reflected pulse has a radius of 5.5 km which 

constitutes a half angle of five milliradians. Consequently, the angle of 

incidence is very nearly zero during this portion of the pulse, with a 

total variation of + 5 milliradians. With the possible exception of a 

completely flat sea surface, this variation of incidence angle can be 

neglected and the analysis can be confined to the case of normal incidence. 

The significance of this simplification is apparent when one considers the 

complexity of an analysis which attempts to predict angular dependence of 

scattered power over the range of angles from 0-90 degrees. Barrack’'" has 

demonstrated that the usual approximations made in evaluating the vector 

Kirchoff integrals for rough surfaces generally lead to results which may be 

in error at angles of incidence greater than 20 degrees. There is, however, 

1 2 

excellent unanimity for the theory near normal incidence; Barrick , Hagfors , 
■ 3 

and Fung and -Moore have shown that models employing geometrical and physical 
optics can lead to identical answers near normal incidence. 

The details of radar backscattering from the sea surface are presented 
in paragraph 4.3, where allowance is made for the fact that the energy 
spectrum of a sea is in general anisotropic. The present knowledge of 
the sea surface energy spectrum is not sufficint to allow calculation of 
these effects. At normal incidence the radar cross-section is known to be 
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a function of the mean-square slope of the sea surface. A few measurements 

of the mean-square slope and the mean-square height of ocean waves have been 
5 

made. Munk observed that when wave numbers corresponding to spatial 

wavelengths of less than one foot were removed from the ocean spectrum by 

means of artificial slicks, the mean-square height of the waves (proportional 

to the integral of S(k)) was essentially unaffected but that the measured 

2 

mean-square slope (proportional to the integral of k S(k)) was reduced by a 
factor of three. This means that although the scattered power at normal 
incidence is heavily influenced by sea slopes, it is rather insensitive to 
the mean-square height of the waves, the latter being a commonly observed 
oceanographic variable. 

4.2 SUMMARY 

With X-band satellite altimetry, for waves with rms heights greater than 
one-half foot, 

1. The backscatter cross-section at normal incidence will in 
general be inversely proportional to the mean-square slope 
and the degree of anisotropy of the sea surface; 

2. assuming a Gaussian shaped autocorrelation function, the 
incremental signals, sampled at nanosecond intervals for 
times up to 150 nanoseconds (after the leading edge of the 
pulse hits the sea) , will be independent; and 

3. the incremental areas uncovered at nanosecond intervals are 
large enough to expect that the distribution of the amplitudes 
of the incremental signals will be Rayleigh-like. 

The significance of these results to the simulation of sea returns was 
discussed in Section 3. In regard to the geodetic altimeter program, 

g 

this study has led to the following conclusions : Backscattered power 

can be. related to mean-square slope of the sea surface. Experimental 
programs using radar scatterometers have indicated that empirical 
relationships can be used to relate the scattered power (slope dependent) 
to ocean surface winds. A considerable extension of the presently 
available oceanographic information (e.g., two-dimensional slope spectra) 
would be needed to place these empirical relationships on a theoretical 
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basis. Assuming that the empirical results can be used reliably to infer 
surface winds, the wind information may be used to infer rms wave height 
forecast programs and other uses. As discussed in Section 5, the oceano- 
graphic relationships between mean-square slope and wind velocity are not 
known to the extent required to permit an electromagnetic solution to the 
problem. The ocean measurement problems are immense. 

4 . 3 RADAR BACKSCATTERING FROM THE SEA 

The design of a satellite-borne radar altimeter will be dependent 
upon knowledge of the characteristics of the radar backscattering from the 
sea surface. Backscattering at all incidence angles from a rough surface 
has received a great deal of attention by numerous authors. The subject 
is most complex and there is no uniform agreement except near normal 
incidence for a very rough surface. 

1 

At near normal incidence, Bar rick has shown that the radar back- 
scattering cross-section is given by 


a = Kk^ exp {-4k^h^[l - p(x,x', y, y 1 )]) dxdy dx'dy 1 


( 1 ) 


where: the limits on the integrals are defined by the area illuminated by 

the radar pulse; 

K is a function of sea surface conditions (mean-squared slope) and 
of the Fresnel reflection coefficient (near normal incidence K is 
relatively insensitive to sea surface slope) ; 

2 

h Q is mean-squared height of the waves; 

p (x, x', y, y T ) is the normalized surface spatial correlation 
coefficient; 

2ir 

k is the Fourier wave number = — in which A is r-f wavelength; and 

(x,y) and (x',y') are positions on the mean sea surface. 

The derivation of this equation is based upon physical optics 
approximations to the Stratton-Chu integral equation. The major assumption 
is the use of the tangent plane approximation to the surface which requires 
that the radius of curvature of the surface be large in comparison to the 
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wavelength of the incident radiation. At a wavelength of three centimeters 
(X-band) it is not clear that the tangent plane approximation is appropriate 
for a sea surface with capillary waves. On the other hand, the capillary 
waves are only millimeters high and, at normal incidence, the effect of 
millimeter scale irregularities on three centimeter energy is not expected 
to be large. It is recognized, however, that there may be some error 
involved in the use of the tangent plane approximation. 

A second approximation has been introduced by assuming that the joint 

probability density distribution of the wave heights at positions (x,y) and 

(x 1 ,y 1 ) is a normal (Gaussian) distribution with zero mean value and variance 
2 4 

h Q . Kinsman , among others, presents experimental evidence which indicates 
that the height of the sea surface is, only to a first order, a normally 
distributed variable. Since the surface statistics were assumed to be 
homogeneous, the correlation of the height fluctuations at any two points 
is a function only of the difference of the coordinates of the points and 
not of the absolute value of the positions of the two points. 

The homogeneous assumption allows (1) to be expressed as follows: 


a = Kk 2 | || exp {- ^^ 2 ' ^ [1 - p (x-x ' , y-y ' ) ] } dxdy dx'dy' 


with x-x* = Ax 
y-y' = Ay 

and integrating oven an L^, square area gives 


L L x+L y+L 
x y x y 

r j£ — r — r -*r- 


= Kk 


2 2 


2 r 16ir 2 h 2 r 1 
exp { x — [1 


- p (Ax, Ay) ]} dAxdAy dxdy. (2a) 


-L -L x-L y-L 
x £ x J 

2 2 2 2 
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1 6 tl 

Note that since — ?>> 1> if P (Ax, Ay) departs from one, l,e,, 

for Ax and Ay differing much from zero, the exponent in the integrand 
would become quite large. The major contribution to the relative coordinate 
integral ( in 2a) occurs for Ax,Ay~0. The exact limits of the relative 
integral are hot important (so long as L and L are large in comparison 

x y 

to the regions of Ax and Ay which contribute to the integral) and we can 
extend the limits of the relative coordinate integral to infinity without 
incurring appreciable error. 


a 



-L -L — < » —05 
_2L_ _JL 


2 2 


exp {- 


2 2 

16ir h r , , „ . \ * 

5 [1-p (Ax, Ay) ]} 

A 


dAxdAy dxdy 


( 2 b) 


2 . 2 


a = Kk 2 Area jj exp {- 16 ^ [1- ( x, y) ]} dAxdAy. 


( 2 c) 


Before proceeding,, it should be noted that a sea surface composed of 
infinitely long-crested waves of a single frequency will cause 
the above arguments to fail. Suppose the waves were traveling 
in the x direction. The correlation function in the x direction 
will be periodic and the relative coordinates in the x direction cannot 
be replaced with infinite limits. In the y direction, the correlation 
function will be unity for all y. The relative limits for y cannot be 
replaced with infinite limits either. While it is true that on a real 
ocean the likelihood of a single frequency, infinitely long-crested 
wave is zero (the ocean is of finite size) , with a sufficiently small 
antenna footprint (corresponding to an extremely short pulse length) a 
moderately long-crested narrow band swell could seriously invalidate 
equations ( 2 b) and ( 2 c) . 

If the correlation function is assumed differentiable at the origin, 
then a Taylor series expansion about the point (0,0) is 
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p(Ax,Ay) = 


p(0,0) + Ax I 
3Ax- j 

Ax 


+ Ay 
a Ay=0 


la 

3 Ay 


IML 


a 2 

Lx 


AxyAy^O 


: Ax: 


Ax=Ay=0 


(Ay) 2 

2 


iffi. 

3 Ay 2 


+CAx) (Ay) a ^ y 


+. 


Ax= Ay-0 


Ax=Ay=0 


(3a) 


Xn the following development we shall show that 


Bp . _ . 3p 

3Ax 3 Ay 


=0, 


Ax=Ay=0 Ax=Ay=0 


consequently. 


where 


Thus, 


p (Ax, Ay) f l-a 2 (Ax) 2 -b 2 (Ay) 2 _ c 2 (Ax)(Ay). 


a 2 = -1/2 -2-2-, 


3Ax 


2-1 

Ax=Ay=0 


b = -i/2 

3Ay 


Ax=Ay=0 


c 2 . -is. 

3Ax 3Ay 

Ax==Ay=0 


(3b) 


1-p (Ax, Ay) 


a 2 (Ax) 2 + b 2 


? 9 

(Ay) + c 


AxAy 


( 4 ) 
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Equation (4) shows that the isocorrelation contours for small values of 
x and y are elliptical- For larger values of x and y more terms 
of the Taylor series expansion would he required and the contours need 
not be elliptical. However, for very high, frequency backscattering, the 
elliptical approximation appears reasonable. The proof of the preceeding 
discussion follows. 

Let z(x,y,t) represent the height of the sea surface at the position 
x,y and the time t 


z(x,y) 


S(k x ,k y ) exp {i[kx+k y y^ dk^ 


From transform theory, since z(x,y,t) is a real variable, S(k x ,k y ) 

S* (-k -k ) . If the height fluctuations are averaged at two points x,y and 

x y 

x' , y' we obtain 


00 

z(x,y) z (x',y ')== 'j 


S (k x »k y ) S (k^, ,k y , ) 


exp r i[k x-k'x'+k y-k'y’l \ 
^ t x x y y / 

dk dk dk’dk’ 
x y x y 


(5) 


For homogeneous, stationary sea surface statistics, z(x,y) z(x'y’) should 
be independent of time t and a function only of x 1 — x and y ,— y« This can 
only be true if 


S(k x ,k y ) S*(k’,k’) - |s(k x ,k y )l 6(k x 




6(k 


■9 


( 6 ) 


where 


isovyi 2 = i s <-yyf 2 


F(k ,k ) ^ F(-k 


x 


x 
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since 


z(x,y) z (x ' ,y ’ ) = R(Ax,Ay) 


fff F(k ,k ) 5(k - k’) 6(k - k*) 

jJJ x’ y x x y y 


exp {i[(k - k’) x - k’Ax + (k - k') y - k’Ay]} 

Y x x x y y J y 


dk dk dk'dk’ 
x y x y 


R(Ax, Ay) = 


F(k ,k ) exp {-i[k Ax + k Ay]} dk dk dk’dk* 
x’y x y J xyxy 


p(Ax, Ay) = 


R( x, y) = Li 

R(0,0) 


F(k ,k ) cos [k Ax + k Ay] dk dk 
x y x y xy 


F(k ,k ) dk dk 
x’ y x y 


Therefore 


9p 
8 Ax 


k F(k .k ) sin [k Ax + k Ay] dk dk 
x xy x y xy 


F(k ,k ) dk dk 
jj x’ y x y 


9p 

3 Ay 


[ k F(k ,k ) sin [k Ax + k Ay] dk dk 
jx x y x y xy 


rr 


F(k ,k ) dk dk 
x y x y 


s z p 


-w 


k k F(k k ) cos (k Ax + k Ay) dk dk 
xy xy x y xy 


3Ax3Ay 


II 


F(k ,1c ) dk dk 
x y x y 


(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 
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a 2 P 

e 

3 Ax' 


f( 2 

k F(k k ) cos (k Ax+k Ay) dk dk 
x x y x y ; i y 


F(k k ) dk dk 
x y x y 


9 Ay 2 


k F(k k ) cos’fk Ax+k Ay] dk dk 
6 y «-y 1 x y x y 


11 


F(k'k ) dk'dk 
x y x y 


Evaluating (8), (9), (10), (11), (12) , (13) for Ax = Ay=0 gives; 


P (0,0) = 1 


Ael 

3Ax 1 
Ax=Ay 


Ifi. 


= 0 


9Ay ] 
Ax=Ay=0 


-1/2 3 2 p 


- -2 
9 Ax 


1/2 


oo 

11 


k F(k k ) dk dk 
x v x y x y 


Ax=Ay 


11 


F(k k ) dk dk 
v x y x y 


= a 


-1/2 9 2 p 


3 Ay 


1/2 


11 


k F(k k ) dk dk 
y x y' x y 


Ax=Ay=0 


F(k k ) dk dk 
x y x y 


= Id 


3 2 P 


3Ax9Ay 


Ax=Ay=0 


11 


k k F(k k ) dk dk 
x y x y' x y 


F(k k ) dk dk 
v x y x y 


= c 


( 12 ) 


(13) 


(14) 


(15) 


(16) 
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It is noted that in view of (8) , the correlation function can be 
considered to result from waves travelling in the k^jk directions. 

Since F(k x ,k^) = F(-k x ,-k ), the correlation function could equally well 
arise from waves travelling in the -k^, -k directions. There is an inherent 
180° ambiguity in the direction of travel of the waves which make up the 
spectrum. 

4.4 CALCULATION OF THE NORMAL INCIDENCE RADAR BACKSCATTERING CROSS-SECTION 


The physical optics approximation to the normal incidence backs cattering 
cross-section per unit area of illuminated rough sea leads to 


co ^ 

a /area = Kk 2 j exp{-(^-?p- ) [a 2 (Ax) 2 + b 2 (Ay) 2 + c 2 (Ax)(Ay)] 


(17) 


dAxdAy 


It may not prove feasible to use oceanographic spectra to evaluate 
2 2 2 

the constants a , h and c since oceanographic spectra do not usually 
have sufficient detail ( at the high wave number end of the spectrum) to 
permit an accurate evaluation through (11), (12), or (13). Nonetheless, 
it is possible to relate the cross-section to an oceanographic variable 
(even though there may be some doubt as to the pertinence of the oceanographic 
measurement) . 

The maximum slope of the sea surface at the point z(x,y) is given by 


Vz (x,y) 


3z ~ , 3z 

S x+ s7 


y 


The mean square slope of the surface is 




3z 

2 

3z 

Vz (x,y) 

2 = 

3x 

f 

3y 


It can be shorn that if z(x,y) has a power spectrum F(k x ,k^) then 


3z 

2 

3x 



and 


M 

1 By 1 


- 2 2 
2 have power spectral densities k F(k ,k ) and k F(k ,k ) 

x x 5 y y x* y 
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respectively. The mean squared values of 3z_ and _3z_ are given by 

3x 3y 


az 

3x 


2 

k 2 F(k k ) dk dk = -h 2 - — 
x ' x y x y 2 

J J 3Ax 


2 2 
= 2h a 

Ax=Ay=0 


2 

s = 

y 

Ax=Ay=0 



ft 


k 


F(k k ) 
x y 


dk dk = 
x y 


3JL 2 

3 Ay 


= 2h 2 b 2 


( 18 ) 


(19) 


Substituting 
a / area = Kk 2 


(18) and (19) into (20) gives 


. 


exp{ 



s 2 s 2 

[ -J (Ax) 2 + CAy) 2 + h 2 c 2 AxAy ]} 


dAxdAy 

2 

The term c can be related to the slope covariance, i.e. 


( 20 ) 



o 


2 

c 


3z 3z | 
3x 3y | 


In general then, the normal incidence radar backscattering is an inverse 

function of the mean square slope of the sea. It should also be observed 

that the cross-section is also a function of the degree of anisotropy 

2 2 

of the sea surface as is evidenced by the term h Q c (Ax) (Ay) . 

One-dimensional ocean wave spectra have been obtained from measurements 
of wave height versus time at a single point . These measurements are 
first processed to yield a frequency spectrum which, through the application 
of the wave dispersion relation, is then interpreted as one-dimensional 
spatial spectrum. This type of analysis is not without its difficulties 
because of the averaging of the high frequency components by the transducer 
as well as the neglect of nonlinearities in the wave equation. 

Two-dimensional spectra can be obtained directly from stereo photographs 
of the sea as demonstrated by Cote, Pierson, et al. 2 . For a variety 
of reasons, not the least of which is the sheer magnitude of the task, the 
spectra so derived do not represent the components whose spatial periodicity 
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is less than 30 feet. Consequently the high wave number end of the 
ocean wave spectrum is not well known. It should be noted that the very 
high wave number end of the spectrum is far less important to the 
oceanographer than it is to the radar user. 

For these reasons the importance of the slope covariance on radar 
cross section has not as yet been evaluated. Frequently, in the absence 
of knowledge about the detailed behavior of the spectral density function 
F(k ,k ) at high frequencies, the radar baclcscattering cross-section is 

“ y 

specialized to the case of an isotropic sea surface 

2 2 
s is equal to s 
x ^ y 

surface 


and c , the covariance 


8z 


9x 

dy 


Under these circumstances 
is zero. For an isotropic 


o I area = Kk 


UJ 

II 


exp { 
dAxdAy 


-8V 


.,2 


le* (Ax ) 2 + s 2 (Ay) 2 ]} 
x y 


( 21 ) 


, Kit 

a / area = — - v, 

s 


Equation (21) is particularly simple and clearly shows the inverse 
dependence of radar cross-section upon the mean-square slope. 

In order to make clear the physical basis for some of the arguments 
which follow, (2c) will be re-examined for the case of an isotropic sea 
surface. In this case 


Cf= Kk 2 area jj exp {— 6 - -j- [i - p(Ax,Ay)]> dAxdAy 


( 22 ) 


2 2 

o/area= 2 Kk 2 fexp{ — — »■ [l-p(r)]} rdr 


where p(r ) 2 = (Ax ) 2 + (Ay) 


= 4 


and p(r) = p W(Ax ) 2 + (Ay ) 2 


which is seen to be the isotropic case. 
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It is clear that if y is a large number, the major contribution 

to (22) will occur for small values of p (r) such that p(r) = 1—6 

whered is a very small number. A two term Taylor series expansion for p (r) 

is given by 


P (r) = 1 


2 a 2 
z ~ 2 


+ . * • 
r=0 


cr /area = Kk^ 2ir |exp ^ (1/2 


.9 ^p I )r^} rdr 
„ 2 I (23) 

or r=o 


cr /area = 2irKk 


2 A 


2 . 2 . 




16tt h (1/2) 3r 2j 
r=0 


(24) 


Comparison of (21) with (24) will show that 
2 




must 


equal Direct calculation will verify that a /area = — . Therefore the 

s 

radar cross-section for an isotropic rough surface is inversely proportional 
to the mean squared slope. 

The significant feature of this last derivation, which will be used 
in certain following sections, is that for a rough isotropic surface 
the major contribution to the integral in '(22) occurs for relatively 
small values of r such that the p (r) - 1. 


4.5 COHERENCE OF BACKSCATTERED SIGNAL 

If the entire illuminated surface of the sea were to act as a flat 
coherent reflector, the returned power should be proportional to the 
square of the area illuminated. However, the arguments following equation (2) 
indicate that when the illuminated area exceeds a critical size, the return 
is proportional to the first power of the illuminated area. Just what 
are the dimensions over which the sea acts as a coherent reflector? The 
linear dimensions can never exceed (at high frequencies) the surface 
correlation length L, and in fact will be shown to be considerably smaller. 
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Consider, for illustrative purposes only, that the normalized correlation 
function of the sea surface is Gaussian: 

2 2 

p(f) = exp {- = 1- f ° r small r. (25) 

L L 

At x-band, with an rms wave height of one-half foot, the analysis 
shows that the main contribution to the integral in (2) occured for small 
values of rfLr^ where r^ is defined by 

[4|k ]2 [ 1- p( Ti )] = 4 (26) 

Therefore, 

p( ri ) - .999 =1- (J) 2 


and 

~ = 3 x 10" 2 

This value of r, which is less than 3 percent of the correlation length, 
is a measure of the distance over which fields are correlated on the sea 
surface. Since the contribution to the integral is quite small outside 
°f r ^ is considered to be a conservative estimate of the linear size 
of a coherent scattering element. Taking r^ to be 3 percent of the surface 
correlation length, for L as large as 500 meters, r^ is still only 15 meters. 

The number of independent scattering areas can be estimated as follows: 
The radar footprint on the sea— surface expands radially with time. During 
each one nanosecond time increment that the radius of the. footprint 
increases by 15 meters, the scattering from the incremental areas shall 
be considered to be independent (uncorrelated) of the previously uncovered 
area. The radius of the radar footprint as a function of time, t, after 
the leading edge of the pulse has reached the sea-surface is given by 

r L = [act] 1/2 for t< t (27) 


where r is the footprint 

a is the satellite altitude 
t is the pulse length. 
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If a is taken to be 1000 kilometers, then for t = 100 nanoseconds 

P . 

the increment of the radius uncovered between (t = 99 to 100 nanoseconds) 


Ar 100 " r 100 - r 99 ' 25 Beters 


for 4r 100 ” r i 


In a similar fashion, = r ^Q “ r i49 a PP rox i mate ly 20 meters. 

Therefore, incremental signals at nanosecond intervals are incoherent 
for times as long as 150 nanoseconds after the leading edge of the pulse 
reaches the sea and the footprint radius is approximately 550 meters. 

4.6 FREQUENCY DEPENDENCY 


The results obtained in this section can be related to the geodetic 
altimeter as follows : An operating wavelength of the satellite altimeter 

at X-band (3 cm wavelength) is currently being considered. Kinsman^ estimates 
that 80 percent of ocean waves have heights greater than three feet. It 
therefore seems reasonable to assume that h, the rms wave height, is 

greater than one-half foot for an appreciable fraction of the time. At 

2 2 2 3 

X-band the quantity 16tt h /A has an approximate value of 4 x 10 . Most 

of the contribution to the integral of equation (22) will occur over a 
range of r from zero tc r^, where r^ is defined by equation (26) . Thus 
p(r) is found to have a value of 0.999. Over this range it appears reason- 
able to approximate p(r) by equation (25) since the change in p(r) is small 
from r = 0 to r^. On the other hand, if the operating wavelength were 
75 cm (400 MHz) the term in equation (26) 


( 


4irh 
o 

A 


6.4 


with the result that the correlation function for the same range of r 
would decay to a much smaller value (p(r^) = .375), 

Therefore, at 400 MHz one would not expect that p (r) could be well 
represented by a simple quadratic. It is concluded that the backscattered 
power would not exhibit a simple slope dependence at 400 MHz if the wave 
heights were only one-half foot rms. At X-band, however, for- the same 
wave, the backscattering cross-section at normal incidence can be accurately 
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related to Inverse mean-square slope,. 


4.7 ANALYSIS OF WAVE HEIGHT EFFECTS ON THE ALTIMETER WAVEFORM 


The extent to which mean-sea-level measurements derived from an 
electromagnetic signal depart from true mean-sea-level is presently 
unknown. If this type of bias is not a significant error source in the 
altimeter data, then estimation of sea state characteristics can be 
considered as an elective requirement in the altimeter program. Other- 
wise, wave height estimation is essential in the satellite. The problem 
of bias in the mean-sea-level measurement is discussed below. 

An extension of the work of Rice and Longuet-Higgins has been formulated. 
In the one-dimensional case, the distribution of reflecting points F(z) 
on a rough surface of height variable z is 


F(z)dz = — 
n 


p(z,z ’ ,z") j z" [dzdz" 


where p is the joint probability of the density of the function z'and its 
derivatives (primed) and n is the normalization term. In order to specify 
that a point on a rough surface constitutes a normal incidence reflection 
zone, the first derivative of the height variable z ’ must be equal to 
zero (e.g. a critical point) and the second derivative z" must be within 
the required curvature range. 

In order to proceed with the mathematical development, it is necessary 
to assume Gaussian statistics. The approach is discussed in Appendix D. 

The degree to which the Gaussian assumption holds is unknown. Figure 4.1 
shows measured probability distributions of the ocean surface for the 
height variable and its slope z*. Note that the slope distribution departs 
much more from a Gaussian statistic than does the height variable. It is 
in fact possible for the height distribution to be symmetrical and for the 
slope and curvature distributions to be skewed. In this case, the radar- 
measured mean-sea-level would certainly depart from its true value. 

It is difficult to pursue this theory because of the above-mentioned 
lack of oceanographic data. The principal oceanographic unknowns involve 
the high frequency behavior of the slope and curvature spectra. With data 
presently available, convergence of the spectrum interval is not even 
assured. 
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Crosswind (upper) and upwind (lower) 
components of slope. Solid curve- 
measured, dashed curve- gaussian. 
From Cox and Munk 


Fig. 4.1. Comparison of measured height and 


0 12 

Units of height squared (z 4 0 ) 
Telative to msl 


lope distributions with gaussian curves 






Because of the assumptions required in a mathematical formulation, 
it is desirable to examine the bias problem using actual ocean height 
profiles in a numerical evaluation of the integration. Any errors that 
may be introduced by the simpler models and assumptions could then be 
evaluated. Such data can apparently be obtained from laser profilometer 
measurements and this information would permit completion of the wave 
model task. Studies could then be conducted to determine (1) feasibility 
of wave height extractions in the leading portion of the altimeter signal 
and (2) importance of "sea-state bias". 
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Section 4 
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SECTION 5 


OCEANOGRAPHIC STUDY 


5.1 SUMMARY AND CONCLUSIONS 

In the course of investigating the effects of mean T^ave-height and 
spectral anisotropy on the scattering of electromagnetic waves, a number 
of oceanographic information deficiencies were encountered. The work 
reported in this section was undertaken as a means of providing statistical 
characteristics of the ocean surface for use in the work reported in the 
previous section. This section, therefore, contains a review of prior 
oceanographic investigations in regard to the dependency of surface 
characteristics on oceanographic and meteorological parameters. 

.The general conclusions reached in this section are as follows: 

(1) The ocean wind waves are probably generated by the 
resonance mechanism proposed by Phillips, but their subsequent 
growth is supported by Miles’ shear flow model. In addition 
to the wind, there are other modification mechanisms, such as 
wave-wave, and wave- turbulence interactions, etc., that will 
change the form and pattern of the wind waves. 

(2) Because of the favored direction and the modification 
mechanisms in the wave motion, the statistical distribution of 
the surface displacement and surface slope are significantly 
skewed . 

(3) By the predominant wind force on waves, the main 
energy containing part in the two-dimensional wave spectrum is 
also skewed. In other words, the wave field itself might very 
well be homogeneous but definitely anisotropic. The 
anisotropy might be gradually smoothed by the modification 
mechanisms, especially at high frequencies, and become small 
but it would never be completely eliminated. 
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Conclusions specifically relating to the integrated one-dimensional 
spectrum are as follows: 

(1) For a well-developed sea, the spectral form of the 
high frequency range changes inversely as the fifth power of 
frequency, or inversely as the fourth power of wave number. 

(2) The spectral maximum occurs at a frequency n Q 
which is equal to the gravity term divided by wind velocity. 

For most practical purposes, n Q could be used as the cutoff 
frequency at the lower end of the equilibrium range. 

— 2 

(3) The variance of surface slope, (V£) is logarithmically 

— 2 

related to surface wind speed. The dependence of (V?) on fetch 
and duration is critical only for a short distance or time 
period respectively. 

(4) With moderate to strong winds (< 40 knots) , the 

— 2 

relation between (V?) and wind speed offers a good method of 
determining one if the other is known. For wind speeds of 
40 knots or less, their relationship is reasonably well 
established. At higher wind speeds the relationship has not 
been verified. 

5.2 BACKGROUND 

The relationship between wind and waves has been obvious, even to 

casual observers, since time immemorial. Early field and laboratory 

studies provided generally unrelated information and hypotheses. The 

fact that the randomness and irregularity of ocean wind waves defied 

description by the earlier investigators led Rayleigh to write: "The 

basic law of the Seaway is the apparent lack of any law." With the 

turning of the century, and especially within the last thirty years, 

statistical theory began to find its way into applications^. The 

2 3 

publications of Longuet-Higgins and Pierson began a new era in ocean 
wind wave study. 
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From the statistical theory, oceanographers learned that in the study 
of random wind waves only the various statistical measures of the motion 
can be regarded as significant observationally, or as predictable 
theoretically. The fundamental measure of an unknown random wave field 
is, of course, the joint probability density of the variables. .The 
second moment of this probability density is of central importance 
because the Fourier transform of it gives the wave energy spectrum which 
measures the energy contained in a particular wave field. The wave energy 
spectrum is of prime importance since it is readily measurable. Naturally, 
the study of the wave energy spectra has become one of the central 
problems of oceanography. 

In the ocean, most, if not all, of the energy of the waves comes from 

wind. However, the detailed physics involving the energy transfer between 

wind and waves is still not completely known, in spite of the important 

works of Phillips^ ^ and Miles^ Numerous observations and theoretical 

studies have revealed some definite statistical properties which represent 

at least part of the whole story, and valuable working relationships between 

some of the salient parameters involved. In order to serve the immediate 

purpose, the properties of the high frequency end of the spectrum will be 

examined first. Although the high frequency waves contain a very limited 

amount of energy, sea surface slope is strongly dependent upon them. 

Knowledge of the surface slope-wind relationship could provide the basis 

11 12 13-14 

for new measurement techniques (Hulburt, Schooley, Cox and Munk 

1C 14 

Cox, and Stilwell ). 

The following is a review of what is known about the statistical 
properties at the high frequency end of the wave spectrum (known as the 
equilibrium range), and one aspect of their application. 

5 .3 THE EQUILIBRIUM RANGE OF THE SPECTRUM 

The notion of an equilibrium range of a spectrum was first advanced 
17—18 

by Kolmogoroff in the study of turbulence. A similar idea was applied 

19 

in an ocean wave study by Phillips . The basic concept of an equilibrium 
range is the same both in turbulence and ocean waves, although the physical 
processes involved are quite different. 
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It is observed that as a turbulent wind blows across an ocean surface, 

energy is continuously transferred to the sea, resulting in wave development,. 

If the duration of the wind is sufficient, some components of the waves 

presumably will continue to grow, while others, especially the high 

frequency ones, will reach a state of saturation — limited by instability 

in the form of breaking. The limiting configuration of a progressive wave 

20 

was studied by Stokes . He found that when the crest angle reached 120°, 

the acceleration of a fluid particle at the crest was equal to one hali; 

the acceleration of gravity, and at this point the crest, or top of the 

21 

wave, would break off. Subsequent observations by Gaillard confirmed 

Stokes' calculation. Ocean waves, however, are irregular with components 

in many directions and of different frequencies and wave lengths. Breaking 

may occur: 1) under the same conditions as described by Stokes; 2) when 

two crests of the wave pattern run together; 3) when a wave moves into 

a region of high energy density; or 4) when short waves riding over the 

crest of longer ones have an energy excess as a result of the radiation 
22-23 

stress . At any rate, it might be reasonably asserted that the wave 

breaks at some locality when the fluid particles have an acceleration 

comparable to the acceleration of gravity (g). Energy is dissipated during 

breaking and at the same time supplied continuously from the wind. At 

some point, the amount of energy supplied just equals the amount dissipated 

19 

and the equilibrium state is reached. Under this assumption, Phillips , 
using dimensional analysis, defined the functional form of the frequency 
spectrum for the equilibrium range: 


<j>(n) = Pg 2 n 5 ; 


n >> n >> n 
Y o 


( 1 ) 


where $ (n) is the frequency spectrum, g is a constant, n 
cps, n^ is the frequency at the spectral maximum in cps, 
is the typical frequency for capillary waves. 


is 


m 


frequency 

and n =(~^— ■)^ / ^ 
Y Y 
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Similarly, the wave number spectrum ip(k) is given as: 


B k 4 

Mk) = (-) ; k » k » k (2) 

— if y o 

where B is a constant containing the angular dependency, k is the wave 

a 1/2 0 

number at the spectral maximum, and k = (“) is the typical wave number 
of capillary waves . 

Various observational results suggest that there is indeed a certain 
range similar to Phillips equilibrium range. Empirical functions fitted 
to the observation indicate a dependence of the -5 power of n or -4 power 
of k. A summary of the dependence of the spectrum on frequency and the 
values of the constant, f3, are given in Table 5-1. 


Table 5-1. Observed Values of Equilibrium Range Constants. 


Observers 

2 — S 

Spectrum n 

Fetch 

Barnett and Wilkerson (1968) 22 

B 

S 



4.95 

0-400 km 

Burling (1959) 28 

1.48xl0” 2 

5.5- 

500-100 m 

Hicks (I960) 29 

1. 21xl0~ 2 

4.9 

16-300 m 

Kinsman (1960) 

1.04xl0“ 2 

4.5 

1,700-3,000 m 

31 

Kitaigorodskii (1962) 

0.65x10“ 2 

5.5 


Pierson (1962) 82 

1.33x10” 2 

-5 

300 km 

33 

Longuet-Higgins , et al. (1963) 

.80x10” 2 

-5 

500 km 


For frequencies or wave numbers greater than n^ or k^ respectively, 

the capillary force would become dominant rather than gravitational force. 

24 25 

The breaking of capillary waves was studied by Crapper and Schooley 

In spite of the different physics involved, the same similarity theory 

26 

still applies. Phillips obtained 


5-5 



( 3 ) 


, v ol 2/3 -7/3 

Kn) = 3'y n ; 


n << n << n 
Y v 


Kk) = (— ) k 

TT 


-4 


k << k << k 

Y V 


(4) 


where 3 ’ and B' are constants, n^ and are the cut-off frequency and wave 

number, respectively, at which viscous forces become important. Observa- 

15 ' 34 

tions in the capillary wave range are scarce, but Cox , ftidy and Plate , 

35 

and Volkov's observations .give some evidence of the Validity of the 
relationship. 

A few words of caution must be added here regarding the use of 
the equilibrium range concept. First, the proposed spectrum function 
represents the magnitude of wave spectrum only for a well-developed sea when 
energy becomes saturated over a wide range of wave numbers, or frequencies, 
manifested by the occurrence of fairly -sharp wave crests and intermittent 
white caps. The sea under this condition is, presumably, statistically 
stationary. During active development or decay, the above relationships 
are not valid. The proposed spectrum, then, can be regarded as representing 
the upper limit of the magnitude of the wave spectrum in these transient 
states. Fortunately, the transient period (or distance) is not very long 
and a quasi-stationary state can be realized within hours with wind 
speeds of 30 to 40 knots or less, or a fetch of approximately 200 km. 

Secondly, the cutoff frequencies of the spectrum are dictated at the 
higher end by capillary forces (in the case of gravity wave spectrum) 
and viscous forces (in case of capillary wave spectrum) . At the lower 
end the cutoff frequency is less definite. However, by considering the 
development of the wind wave spectrum (Phillips^’ ^ ’~^)» the spectrum can be 
taken for most practical purposes as the frequency (or wave number) of the 
spectral maximum. The spectral maximum should occur, according to Phillips' 
resonance theory^ 5 3,26,6 afc t ^ e f re q UenC y or waV e number having the 
following relation 


n o=k . U (5) 
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where n^ and k^ are at spectral maximum, and U is the wind velocity. 
Using the relationship between n and k; i.e. n^ = gk, then 


Hence 


or 


2 2 2 
gk = n « k U . 
° o o o 


k « 
o 



n 

o 


oc 

U 


( 6 ) 


(7) 


This theoretical result was subsequently confirmed by the observations 
of Kitaigorodskii^, Moskowitz^, Pierson^, and Volkov^. Pierson' 
results are summarized in Table 5-2. 


Table 5-2. Some Characteristics' of Pierson's 



(1959) Spectrum. 

A ' 

U cm/sec t hrs. f max. 

A 

U fmax/g 

1 

110 

10.5 

0.11 

0.012 

2 

106 

13.5 

0.11 

0.012 

3 

148 

10.5 

0.08 

0.012 

4 

148 

13.5 

0.08 

0.012 

5 

, 143 

16.5 

0.07 

0.011 

6 

175 

14.5 

0.07 

0.012 

7 

171 

17.5 

0.06 

0.011 


A 

NOTES: 1. U is the frictional velocity; it could be related 

to U, the wind speed, directly. See for example 
Pierson and Moskowitz^l. 

2. t = duration of wind action. 


Besides these observations, there are other attempts to develop an 
empirical expression for the wave spectrum function. Among the better 
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known earlier ones was that due to Neuman , who proposed the following 
expression using wind velocity as a parameter: 


4> (n) = cn 6 exp { — r } (8) 

n IT 

Aside from other objections, the inconsistency of the dimensionality 

makes it unlikely that it is representative of a genuine physical law. 

43 44 28 

Numerous other attempts have been made by Bretschneider ’ , Burling , 

45 46 

and Darbyshire ’ , to cite a few. But at that time the inclusion of 

such a parameter as wind velocity by curve-fitting based on limited data 
seems premature at best. 

No major breakthrough was reported until Kitaigorodskii and 

38 47 38 

Kitaigorodskii and Strekalov ’ published their findings. Kitaigorodskii 

assumed that the energy containing part of the wave spectrum is a function 

of only four variables; i.e. 


Kn) = F(n,g,U*,X) 


(9) 


L T / O 

where U A = frictional velocity = (— ) x/ and X = fetch. 


Then by dimensional analysis, Kitaigorodskii obtained the expression 


U * n 


<f>(n) = g 2 n 5 F 1 (- 


£X_ 


When the expression (- 


U*n 


, -®y) is normalized it does indicate that 
U 


5 41 

«j> (n) n approaches a constant value. Later, Pierson and Moskowitz 


2 

g 


made a slight modification of Kitaigorodskii' s expression by simply 
changing frictional velocity, U A to wind velocity,- U, and arrived at 
the expression: 


( 10 ) 
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( 11 ) 


2 % 4 
$(n) = a -^r e ^ n ^ 
n 

where a , 3 are dimensionless constants, n is the frequency at spectral 
maximum and is equal to -jj-. In both cases, <j>(n) -*■ dg n as n » n^. 

This conclusion gives considerable support to the Phillips equilibrium 
range theory. 

5.4 SLOPE SPECTRUM 

Similar to the energy spectrum discussed in the last section, other 

spectral functions may also exhibit an equilibrium range. Considering 

the frequency spectrum of the surface slope [S R (n)] at a fixed point, as 

50 

an example, Phillips obtains 




2 

7T 


f H(x.t’) 8£(x.t* + t) 


3x 


3x 


cos nt dt 


$ 


( 12 ) 


where x q , Xg are the horizontal vector components and £ is the surface 
elevation. In the equilibrium range, by dimensional analysis. 


S rt (n) = c n n \ n >> n >> n 
a0 ag Y o 


(13) 


where c is a constant tensor. Under the assumption of statistical symmetry 

33 

with respect to the wind direction, c was determined to be 

act 


c = c + c„ — 0.8x10 
aa 11 22 


-2 


(14) 


In the capillary range. 


S „(n) = c T _n 1 , n << n <<n 
a3 aB ’ y v 


( 15 ) 
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and c' is not necessarily equal to C^. An experimental study by Cox^ 
gives some support to the prediction. 

Since the slope spectrum is derived from the product of derivatives 
of the surface elevation, the relationship between them can be shown 
schematically as' in Fig. 5.1. Accordingly, the variance of surface slope is 
given by 


(V?) = I k <Kk; dk 

k 


<j> (n) dn 


n 


(16) 


Equation (16) shows that the contribution to (V?) is heavily weighed 
2 4. 

by the factor k or n in the integrand. Thus a major contribution comes 

from the high frequency (or wave number) end where a negligible amount of 

energy is contained. However, due to the results of recent research in 

optical methods by Schooley"^, Cox and Munk^^’^^, Cox^^, Stilwell^ and in 

27 

radar devices by Barnett and Wilkerson , both using the surface reflections, 
the importance of the study of surface slopes has markedly increased. In 
the following section an important application is presented. 


5.5 THE RELATIONSHIP BETWEEN MEAN SQUARED SLOPE AND SURFACE WIND SPEED 

48 

In the first attempt to optically study ocean surface slopes, Hulburt , 

using sun glitter, concluded that the reflecting facets of the sea surface 

occured most fr equent ly at about 15°. No direct relationship between the 

2 

wind vector and (V?) were found, but it was stated clearly that the wind 

2 

was the major cause of the change of (V?) . 

13 14 48 12 15 

Twenty years elapsed before Cox and Munk ’ ’ Schooley , and Cox 

conducted another series of experiments. The principle used in all the 
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RELATIVE SCALE 



WAVE NUMBER K 


Fig. 5.1. Sketch of Elevation (Energy) and Slope Spectral Distributions 



optical measurements of the surface slope was the same, but t he impo rtance 

• 2 

of the latter studies lay in a proposed relationship between (V£) 

and the surface wind. 

12 

Schooley used flash photography taken from a bridge. He made no 

serious attempt to find an empirical formula to describe the observations. 

At about the same time period, Cox and Munk 33 * ^ studied sun glitter 

photos taken from an airplane. The cases of a clean surface and a surface 

covered with artificially induced oil slicks were bo th studied. Their 

2 

extensive measurements led them to claim that (V£) varied linearly with 

wind velocity. In an attempt to lend more support to their curve fitting 

result, Neuman's spectrum was used to calculate the relatio nship approximately. 

2 

From the calculation, a linear relationship between (V?) and U was indeed 

—6 

obtained. However, since Neuman erroneously used ^(nj^n for high 
frequency, the value of the calculation and the result are questionable. 

Later Cox 33 made a laboratory study and concluded that the relationship 
was not a simple linear one. Furthermore, the value becomes independent 
of fetch after a short critical distance. These findings are in general 
agreement with Phillips' equilibrium range theory. 

Since Phillips’ equilibrium range theory is consistent with numer ous 

2 

independent observations it was used to derive a relationship between (V?) 

and U. As discussed in Section 4.3, the cutoff wave number could be taken 

as k and k - . Then 

c o u 2 




k (k) dk 


k 3 (k) dk = B log (^) (17) 

o 


g 1/2 

If k is much less than the capillary wave range; i.e. k << k = (■“■) , 

c , c y y 


then 


k IT 
c 


(V?) = B log (-— — ) . 

O 


( 18 ) 
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then 


However , 


if k 


c 


is within or greater than the capillary wave range 


(VC) 2 « B* log (^-) + B log (^) (19) 

Y ° 

When these calculations are plotted together with data from Cox 15 and 
13 14 48 

Cox and Munk * * for the case of an oil slick which filtered out 

almost completely waves with length less than 30 cm- the Phillips 

calculation shows remarkable agreement. For the case with capillary 

components, the agreement becomes less as wind v elocity increases. This 

2 

could be due to the difficulty in measurement of (VC) as reported by 
Cox and Munk . 

41 g 

Using Pierson! s spectrum and integrating from n^ = ^ to some 

frequency n, with the help of the dispersive relationship, gk - n, 'then 


(V?) Z 



2 4 

5 2 
m g 


n > 


dm 



o v 

U k 1 , 1 _0 1 

g "2 2 fc 2 


( 20 ) 


This is of the same function form as (18). Since Pierson's formula 

applies only to gravity waves, its usefulness is limited. 

2 

The relation between (V?) and U seems to be well defined both 
experimentally and theoretically; however, the application of this 
relationship is not without difficulty, especially when the wind 
velocity becomes very high (40 knots or higher) . The difficulties are 
twofold. First, as the wind increases, the stability of the wave is 
controlled not only by breaking but also by blow-off from the crest. The 
validity of Phillips 1 equilibrium range theory under this condition has 
not been carefully studied. Secondly, with increased breaking of the 
energy-saturated waves, optical observation becomes more and more 
difficult. White-caps tend to result in an over- estimation of the surface 
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slope . This, however, could be overcome by improving measuring 
techniques. (Clearly, the radar technique is a promising one.) 

5.6 DIRECTIONAL CHARACTERISTICS OF OCEAN WAVE SPECTRA 

Based upon the statistical approach to ocean wave study, the free 
surface is represented as a sum of independent components of sinusoidal 
waves of different amplitudes, phases, and directions of propagation. 

The mathematical ideas involved in this representation are simple enough, 
yet the inverse of the problem; i.e., to find those components which 
accurately represent the random wave field, poses an extremely difficult 
task for oceanographers. This, in essence, is the problem of finding a 
two-dimensional wave spectrum for a random wave field. 

The difficulties encountered in the measurement techniques necessary 
to obtain good two-dimensional spectra remain almost insurmountable. 

However, the practical need for this information has forced researchers 
to seek a modus vivendi ; i.e., to seek a limited goal of determining a 
one- dimensional spectrum to represent an integrated effect of the actual 
two-dimensional phenomenon. This one-dimensional spectrum could be 
obtained by means of a Fourier analysis of a measurement record at a 
single point. Even a one-point measurement was difficult to obtain in 
the deep water a decade or so ago, because of the lack of a working platform. 

The information contained in a one-dimensional spectrum is, indeed, 
valuable as far as the total energy of a given sea state is concerned. 
However, this information is far from complete. The directionality is 
necessary not only from the practical aspects of wave prediction, mass 
transport, impact of waves on coastal structures, etc., but also in the 
basic study of the detailed physical processes involved in wave genera- 
tion. The skewness of the surface distribution has important bearings 
on radar and optical studies of the ocean surface itself. 

A reasonably accurate characterization of these parameters has 
evolved during the last decade through the persistent effort of many 
researchers. These field studies have gradually revealed the statistical 
properties of the ocean waves. The purpose of this discussion is to review 
these efforts, and some of the salient results of the directional statis- 
tical properties. These properties can be divided conveniently into two 
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categories: 1) the probability structure of the ocean surface and 2) the 

directional distribution of the energy in a two-dimensional spectrum. 

The first deals with the geometrical shape of individual waves; the second 
deals with the combined effects of the individual waves upon the wave 
spectrum. 

In order to understand these results, the processes of generation' and 
the mechanisms of modifying the waves will be discussed first. The process 
of generation of waves on the ocean surface are summarized In the next 
section. 

5.7 THEORY OF WAVE GENERATION BY WIND 


The interactions between the atmosphere and the ocean appear in 
various forms. From the oceanic viewpoint and with the atmosphere as the 
primary energy source, motion ranging from global circulations to micro- 
scale turbulence are found. The progress in the study of ocean circulation 

49 . 50 

has been summarized by Stommel and in the study of waves, by Phillips 

51 

Only a few years ago, Ursell , in a review paper on the theory of wave 
generation by the wind, stated that all of the three basic elements for 
determining the mechanism of wave generation were still wanting; namely, 
a correct theory, a well-documented experimental study, and a demonstra- 
tion of a good agreement between the two. This pronouncement served as 

a great stimulation for the study of this problem. Within a year, two 

52 53 

independent theories by Phillips and Miles were advanced. Since then, 

both theoretical and experimental works have resulted in tremendous progress. 

A brief outline of the theoretical aspect is given below. 

52 54 

Phillips * postulated the resonance mechanism of atmospheric 
turbulent pressure and the surface waves. The air flow is always turbulent, 
the pressure fluctuation in the air is, therefore, random. However, as 
the turbulent air flows over a random ocean surface, interaction occurs. 

The fluctuation of the pressure can be separated into two groups. The first 
consists of the purely turbulent pressure fluctuations whose pattern is 
everchanging as the eddies in the turbulence grow, interact and decay. 

The second consists of a spatial pattern of pressure variation that travels 
over the surface at a pressure pattern convective velocity. The existence 
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of this latter pattern was confirmed by observations of Willmarth and 
55 

Wooldridge . To this semi-rigid corrvecting pressure fluctuation, Phillips 
assigned a relationship of the same wave number k as the surface waves. 

The total pressure field thus could be expressed in Fourier-Stieltj es 
representation as 


dm (k,t) = (v + iy) pmc kdA (k,t) + dm(k,t) 


( 21 ) 


where the pressure p(x,t) in the physical space is 

p(x,t) = 

and the surface elevation ?(x*t) is 


. ik*x 
dfo(k,t) e r ~ 


( 22 ) 


£(x-t) 


' dA(k • t) e^‘5 
k 


( 23 ) 


and dA(k*t) is the purely random part uncorrelated to surface waves, 

v,y are the coupling coefficients, c is the phase velocity of the surface 

wave with wave number k, p is water density, and k is |k| . 

to *** ^ 

Through the linearized equations of motion, Phillips was able to 


relate the surface wave spectrum, (j>(a,0) to the pressure spectrum ir(k,x) as 


<j>(a.e) - 


k 2 gt 

2 <Sp/ o 


u 


ir(k,x) cos [ ( -- ■ ■ c r . ° s — 1 ? - l)ax] dx 


( 24 ) 


where x = time variable, t = duration of wind action, U = convective 
velocity of pressure pattern approximately equal to the mean wind velocity. 

There is sufficient evidence to support this theory. First, if there 
is a rigidly convecting pressure pattern, and if the response mechanism is 
indeed responsible for the wave generation, then waves of frequency n 
would be expected to appear in a bimodal distribution in the two-dimensional 
wave spectrum at angles +a dictated by the resonance condition. The 
resonance condition is 


n = k 



cos a 
m 


( 25 ) 
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or 


a 

m 


-1 

cos 


cjj-) = cos 


c-V /2 

kU 


Field observation by SWOP indicated that there was indeed a bimodal 

distribution. Secondly, if resonance did occur, the components of waves 

grow and the fastest should travel at approximately the same speed as 

the mean wind, and the initial growth rate should be linear. This is 

27 

also supported by field observations of Barnett and Wilkerson . 

Complementary to Phillips' resonance theory, but physically quite 

7 

different, is the model suggested by Miles and augmented by Brooke 
58 59 26 

Benjamin , Lighthill and Phillips , Miles considered shear flow over 
existing waves. The mechanism involved is the momentum transfer from 
the perturbed shear flow of the wind to the perturbing waves; therefore, 
this model depends crucially upon the wind profile over the waves. 

Based on this physical picture. Miles calculated the flux of momentum from 
the wind to the waves through a "matched layer", (where U(z ffi ) = c with 
U(s) = mean wind velocity, c = phase velocity of the existing waves) as 


(26) 


T 

0) 


fU 2 , U'\,| 

* p a [ !T ( 'l 


m 


(27) 


where p = air density, U.(z) = mean wind velocity profile, and m 2 = mean 

cl 

square vertical induced fluctuation of the air. 

Since the energy density E of the wave is related to its momentum M by 


E = Me 


the energy flux formula is 

1 - v * Ip « c [ i- 1 1 M • < 23) 

m 

60 26 

Subsequent works of Miles and Phillips uncovered a link between 
the two different theories. From Miles’ model one can calculate the 
coupling coefficients between air pressure fluctuation and wave as 
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V 


(29) 


P a k t 
2 

P c 
u o 


✓ TI v 2 -kz , 

(U cos a - c) e dz 


U = 


2t' (6) 
to 

12 2 
p c k a 
(n 


Since v is definitely negative, the surface pressure should be 

exactly out of phase with the surface elevation. Confirmation by observa- 

33* 61 

tions were reported by Longuet-Higgins et al , and Shemdin and Hsu . 

From the experimental results, it is believed that both Phillips’ 

and Miles’ theories are correct to some extent. Phillips’ resonance 

theory was successful in initiating the wave motion, but it is not an 

effective mechanism to feed in energy to support the continuous growth. 

Miles’ model, on the other hand, depends upon the pre-existing waves to 

perturb the shear flow; therefore the trigger mechanism is missing. 

These short-comings were partly taken care of by subsequent works. Yet, 

another difficulty of the inability to incorporate the fully turbulent 

air flow in the theory is still unsolved (see, for example, Phillips^) . 

In spite of this, Phillips' and Miles’ theories have successfully explained 

many details in wind wave generation. 


5.8 OTHER MECHANISMS OF WAVE GENERATION AND MODIFICATIONS 

Of all the factors having dynamical influence on waves, the primary 
one is, of course, the wind, which is responsible for generation, growth, 
and even under special condition, decaying of the waves. However, the 
wind is by no means the only mechanism that has dynamical consequence. 

Once the energy is fed into the ocean, other processes become dominant 
such as turbulence resonance and shear flow instability. These processes 
cannot create energy, but they certainly transfer energy from one form 
to another. Their results may be shown by a simple change of form of 
waves, and ultimately change the probability structure. Some of these 
mechanisms are discussed below. 

(a) Wave-Wave ' Interaction (nOn-resonant) 

This involves the interaction of waves with their own higher 
harmonics; or in other words, the higher order solutions of the classic 
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wave equation such as Lamb , Pierson and De . The higher order 
contributions tend to modify the geometry of the wave structures, 
(b) Wave-Wave Interaction (non-linear resonant) 

While the bounded interaction discussed above only results 
in the change of waveform, the non-linear resonant interaction 
changes the wave pattern. It transfers energy from an existing 
frequency to previously non-existing ones. The mechanism was 


37 

discovered by Phillips 
and augmented by Phillips*^’ Longuet-Higgins^^, Benny^ 2 , 


for gravity waves, and then simplified 

, I 
70 


68 69 

Hasselmann 5 and Bretherton' w . A similar resonance phenomenon 

71 

in capillary waves was- discovered by McGoldrick 


principle of this resonance is simple, 
waves expressed as 


The mathematical 
Assume systems of interacting 


±n 

e- L 

y=±l 


iX r 
a e : 
r 


X = k 
r ~r 


x - n t 
~ r 


±m k z ix 

4> = 2 be r e r + <f>, (x*z,t) 

Y =±l 


(30) 


I 


where <j>^ is the product of the interaction, ? is free surface eleva- 
tion, and <j> is the velocity potential. Using these expressions in 
the first and second order equation, we can only get similarly 
bounded resonance as in part (a). However, when the third order 
equation is used; i.e. 


ifi 

3t Z 




Mi_ , a_ A + e I c 2 il. 

' " ^ 3z \ 2 S dz J 2 .2 

3t 3z 


3z 


, 3 2 <fr , a ii'. 


3t 


3z 


(31) 


3 2 3 2 

37 u - ? 3in u - 5 


V (y U 2 ) at z = 0 


the interaction of three primary waves of r = 1,2,3 will produce 
a fourth wave of wave number and frequency and respectively 
provide the resonance condition 
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?! + ? 2 + K 3 + K 4 = 0 
\ ± N 2 ± N 3 ± N 4 = 0 


( 32 ) 


Quantitative measurements of resonance interacting waves were 

72 73 

conducted by McGoldrick et al and Longuet-Higgins and Smith 

Their results confirmed Phillips' theory. 

In the case of a random ocean where there are infinitely many 

chances to form the resonant quartets, the interaction mechanism 

will produce wave components of all wave numbers of small energy 

density, including, of course, those traveling against the wind. 

(c) Scattering of Gravity Waves by Turbulence 

Scattering of waves by turbulence has been studied by many 

Ik 

authors. Some of the previous works were summarized by Batchelor 

75 

The true application to oceanic environment was advanced by Phillips 
Physically, it should be fully expected that the presence of 
random velocity fluctuations in the water of a velocity scale 
comparable with particle velocity associated with wave motion and a 
length scale comparable with the wavelength of existing waves will 
result in the convective distortion of the wavefronts, and so the 
establishment of a scattered wave field. In the ocean, it is well 
known that because of wind-induced currents, waves, and wave 
breakings, turbulence is always present in the surface layer. But 
the lack of information on the turbulence on the ocean surface 
layer seriously hampered detailed analysis of the interaction 
between turbulence and waves. As a result, Phillips was forced to 
consider only the weakly turbulent case when 


U t « c(k) (32) 

where U is the r.m.s. turbulent velocity and c(k) is the phase of the wave 
of wave number k. Under this restriction, Phillips found that the 
spectrum of the scattered waves s(k*8).w±ll be 
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s(| • 0) = — '^T Jl ^CK) (33) 

(Sk) i/Z 

where k is the scattered wave number, lc is the incident wave number, <|>(K) 
is the wave number spectrum of turbulence and K is the wave number of 
turbulence . 

Furthermore, the scattering is most effective for a given 0 when 
the scattered wave will have the magnitude of wave number equal to 
the incident wave. 

Unfortunately, due to the difficulty of measurement, no 
well-documented evidence of any sort could be found either to 
support or to dispute this result. However, since the physical 
model used in this analysis is resonable, the analysis is rigorous 
and there is no reason why the result should fail or change 
drastically from that predicted even in a strongly turbulent case. 

Here again a mechanism is present which will produce weak waves 
traveling against the predominant direction of wind force. 

(d) Others 

Aside from the mechanisms discussed, other processes such as 

76 

generation of parasite capillary waves (Longuet-Higgins ) , short 

waves on current or riding on long waves (Longuet-Higgins and Stewart^), 

78 

attenuation of waves by breaking (Phillips ), second-order resonant, 
interaction between capillary waves (McGoldrick^" 1 ") all contribute 
to the change of waveform. But since their influence on the 
directional distribution of energy, or on anisotropic properties 
of the surface, is much smaller than the others, a detailed summary 
is omitted. 

5.9 PROBABILITY STRUCTURE OP THE OCEAN SURFACE ' 


The ocean surface can be regarded as the sum of a large number of 
independent wave components of different wave numbers, phases and 
amplitudes. If the area under consideration is small compared to the 
generating storm area, its statistical properties could be taken as 
stationary . Furthermore, by the central limit theorem, the probability 
distribution of the surface displacement should be Gaussian. Then 
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p(C) = 


(34) 


6X5 i”Cl/2)C 2 /? 2 }. 

(2mS ) 

Field observations invariably indicate that to a first-order 
approximation, the probability distribution is indeed Gaussian as 
reported by Longuet-Higgins ^ 5 • However, because of the favored 

direction of the wind force, and linear and non-linear wave-wave inter- 
action, skewness is always introduced. Therefore, there are some higher 
order modifications of the surface. For example, the third-order 
correction to the Stokes irrotational wave indicates pointedness of the 

crests. Careful field observation by Kinsman 3 ^ indeed substantiated this 

81 

prediction. Longuet-Higgins calculated the correction of the distribu- 
tion from a Gaussian expression and found that the surface elevation 
could be better approximated by the successive terms of a Gram-Charlier 
series as 


p(5) = - i 1/? exp {(1/2) t 2 Hl + (1/6) £ 3 H- + . . .} (35) 

(2 O 1 ' 

where 

2 2—2 3 

t = C /? and H = t - 3t. 

On the other hand, the non-linear interaction and the favored direction 

of wind force also changes the surface slope. Although the non-linear 

37 81 

effect studied by Phillips and Longuet-Higgins showed that the 

—3 _2 3 /2 

change in skewness m^ = C /(£ ) was of the same order as the fourth 

power of surface slope, the wind effect has never been accurately calculated. 

13 48 

Detailed field measurements reported by Cox and Munk * showed the 
skewness much higher than could be accounted for by non-linear effects 
alone. 

Considering the random wave field as a whole, the cumulative effects 
of the wind are manifested clearly in the directional energy distribution. 
Inspection of the two-dimensional energy contours indicates that the 
direction of wave energy distribution is skewed toward the wind direc- 
tion. Even at high wave numbers where there is a considerable amount of 
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scattering and smoothing, the skewness is still obvious. 


5.i0 METHODS OF MEASURING DIRECTIONAL SPECTRA AND SOME RESULTS 

Since oceanographers first adopted a statistical approach to ocean 
wave study, the effort to actually describe a given random ocean by a 
complete spectrum has continued. Unfortunately, the task has never been 
carried out with total success. The difficulty is twofold. Firstly, the 
difficulty arises from the basic statistical and mathematical tools, 
borrowed in whole or in part from other fields. Ideally, by the oceano- 
grapher's definition, a complete wave spectrum of a given sea state will 
describe not only the total energy and the directional distribution of 
energy but, more precisely, the sense of the directional distribution. 
Secondly, by rigorous spectral analysis, under the assumption of 
homogeneity, the sense of the directional distribution of energy cannot 
be determined uniquely. This can be explained as follows. Take the 
surface elevation at position and time t as ^(x • t) , then at position 
x + r, and time t + t, we have (x + r, t + t) . Their correlation function 
R(x,t; r,x) is 


R(x, t; r ,x) ** ?(x,t) ?(x + r, t + x) (36) 

Under the homogeneity assumption, the correlation function R(x,t; r,x) 
should be independent of the origin of the coordinate system, therefore 
R(x,t; r,x) reduces to R(r,x). By a simple translation of coordinates 
in space and time. 


R(r,x) = R(-r, t) = R(r,-x) = R(-r,-x) 


i.e. R(r,x) is an even function with respect to both r and x. 
Mow since the spectrum function i/i(k,n) is defined as 


Kk,n) 


' f Str.T) e iq ?'5- nT) dr dr 
. 

r x 


(37) 


(38) 
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note that ^(k,n) has to be an even function with respect to k and n. 

In other words, there is no difference between 

t^(k,n) and tK-k,-n) or iK~k,n). 

Therefore, strictly speaking, we can get only directional information which 
does not include the sense of the direction from traditional spectral 
analysis. Of course, it must be understood that the ambiguity occurs only 
when the homogeneity assumption is used, but whether this assumption really 
holds is not known exactly. However, without this assumption, or if the 
condition of homogeneity does not exist, then the data collected would be 
representative of that particular place and time. At any rate, the infor- 
mation contained in this conventional spectral analysis is not complete. • 

In order to resolve this dilemma of the it ambiguity under the homogeneity 
assumption, additional physical information has to be added independently 
to the spectral analysis to specify the exact direction and sense which is 
essential to the oceanographer's definition of a complete directional 
spectrum. In simple cases such as in the study of swell or of waves with 
well-defined crest lines, this additional physical information is easy to 
obtain. However, in complicated random wave fields, where individual 
waves are difficult to follow, the true direction of the wave propagation 
can be known only if a complete time history of the whole wave field is 
known. None of the measurement techniques developed to date will provide 
the required information. 

In spite of this difficult problem, some observations have been made, 
limited, of course, by practicality and loaded with assumptions and 
approximations. The following Is a summary of the methods used and results 
obtained. 

(a) Stereo Photography 

A photographic record of the instantaneous state of the sea 
will undoubtedly contain the necessary information needed for deter- 

82 

mining the directional spectrum. This idea was first used by Barber 
to obtain a qualitative estimation of the directional spectrum. 

Later the New York University group embarked on an extensive stereo 
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wave observation project"^ and produced a series of complete 

83 

evaluations of directional spectra (reported by Cote et al and 
84 

U'beroi ). Despite the considerable trouble taken, the results 

obtained were not the real directional spectrum, but instead 

[<j)(a,9) + <j>(a,0 + tt) ] , because of the built-in ambiguity of ir in 

the spectrum analysis. This forced the investigators to assume that 

the energy distribution was confined between (tt/2 ,-tt/ 2) . This 

assumption easily solved the directional ambiguity, and -could be 

accepted as a first approximation for the energy containing range. 

Unfortunately, in the light of the wave modification mechanisms 

previously discussed, a fraction of waves with small, albeit non-zero 

amounts of energy at high wave number can propagate against the mean 

wind. In addition, there even exists the possibility of residual 

low wave number swells coming from other possible directions 

including the one against the wind. Though their omission might 

not be crucial in energy considerations, the high wave number part 

is important in determining the slope and curvature spectrum of the 

sea; therefore, the result obtained by limiting the energy contained 

between (ir/2,-Tr/2) is of questionable quality. 

(b) Directional Array of Probes 

The idea of using a linear array of probes for directional wave 

85 86 

measurement was suggested by Barber * , but the only example of 

measurement was for a single band of frequencies. The same idea was 

87 88 

later picked up by Konyayev and Dreyer and Dreyer and Konyayev , 

but the results so far are still qualitative. 

A slight alteration of this linear array of probes was to arrange 

89 

the probes at the vertex of a triangle. This was tried by Munk et al 

90 

in observing long swells, and by Tsyplukhin in observing near shore 

shallow water waves. In both cases, the waves were relatively simple 

with the fixed direction of propagation; therefore, the results were 

not representative of a truly random wave. 

91 

Krylov et al used the same setup of Tsyplukhin in relatively 
deep water to obtain information on the angular energy spectrum which 
is defined as 
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K6) = 


i|i(a,0) da 


(39) 


where <Ji(a,0) is the two-dimensional spectrum. They found that the 
best fit directional angular spectrum is a cosine power law such as 


<j>(0) a cos (0 - a) 


(40) 


where a is the principal direction of the spectrum. Later, Drylov 
92 

et al proposed an empirical formula for the directional spectrum; 


iKa,6) = 


C 2 t 


2tt 


O r0.32 , 1.8 , 0 . 8 ^ 

i — (cos 0) — — exp ( t— ) } 


(41) 


where t is the average period. Comparison with observation is 

° - 19 

favorable, but compared to Phillips' equilibrium range theory the 

dependence of 1 /t q coupled with the inconsistency of dimensionality 

in. the expression itself would undoubtedly limit its usefulness. 

The idea of using the probe array method is theoretically sound 

but impractical unless stable working platforms and large numbers 

of densely deployed probes become easily accessible in deep water; 

this severely limits its usefulness. 

(c) Buoy Techniques 

82 

This idea was also first suggested by Barber , but subsequently 

2 

developed by Longuet-Higgins and put in practice by an NIO group. The 

instrumentation, measuring process and results were reported by 

80 93 94 

Longuet-Higgins , Longuet-Higgins et al , Cartwright and Smith 

95 

and Erwing . The principle involved using the freely floating buoy 
to measure £, 3?/3x and 3£/3y by gyroscopes, then relate these quantities 
to the first five Fourier components of the spectrum function. In 
this conversion, different combinations of coefficients will produce 
slightly different results, with the same effect as using different 
weighting functions in averaging, but comparisons with the data 
suggested that cosine power laws gave the best fit. The results so 
obtained are the most detailed and complete to date. Some of them are 
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rendered in contour diagrams in Figs. 5.2, 5.3, and 5.4. 

Several points should be brought up here. First, an interesting 
comparison with theoretical prediction on resonance angle could be 
made here. Consider the integral 

1 f 2 2 ®”®2 

I = -7£r J 16 sin ~ 2 — sin “ — 1 K a > 9 ) ^0 (42) 

with ~ arbitrary fixed angles. We can find a set of ^ or 

each given i|»(a,0) to minimize I. It can be shown that 

0 = 1/2C6-L + e 2 ) 

is approximately the main direction of the spectrum, while 

* = i/ 2 ( 0 1 - e 2 ) 

is approximately the r.m.s. angular width of the energy distribution. 
Comparison of with Phillips’ resonance angle shows the experimental 
points are low but not inconsistent with the theoretical result. 

Secondly, there is no clear bimodal distribution of energy, instead 
the best approximation is a cosine power law as 

<f> ( a , 0) “ (cos l/28)^ s (43) 

where s varies with frequency. This power law dependency was 

93 91 92 

reported both by Longuet-Higgins et al and by Krylov et al ’ 
through different methods of measurements. From the form of Longuet- 
Higgins' expression, it can seem that as s decreases so does the 
dependence of <f> on 0 . For a special case of s = 0, the distribution 
of energy will be isotropic. At high frequencies the value of s is 
indeed very small. Therefore, the angular range of energy distribution 
could be expected to be fairly wide. A word of caution should be added. 
Although both Longuet-Higgins and Krylov suggested the cosine power 
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Fig. 5.4. Two-Dimensional Spectrum 


law, and both their power indices decrease with increasing frequency, 

there is a cruicial difference. In Krylov's expression, energy distribu- 

tion is confirmed in a half plane bounded by (y , as implied by the 

dependence of cos (g-a); while in the Longuet-Higgins expression the 

0 

dependency is on cos thus energy could be coming from the other side 
of the half plane at high frequencies. Considering the modification 
mechanism the Longuet-Higgins form seems more realistic. 

Thirdly, the energy spectra contour diagram shows that skewness toward 
the wind directions is apparent, but that energy of small intensity is propa- 
gating against the wind. This upwind propagation of energy could easily be 
accounted for by one or more of the mechanisms discussed earlier. The 
skewness could be attributed to the direction from which the wind is blowing 
and also to the attenuation effect under adverse winds. 

Fourthly, the integrated spectrum over all angular contributions shows 

\ 

the existence of an equilibrium range. 

(d) . Measurements of Orbital Velocity of Waves and Pressure 
and Their Use in Determining the Directional Spectrum 

96 

This idea was first suggested by Nagata and later developed by 
97 98 

Bowden and White and Simpson for measurements in England. The method 
of deriving the directional spectral information from velocity and 
pressure measurements is very similar to the buoy technique. Assuming 
there is a potential function $ for the wave motion. 


u 



3y 


and p = 


1 

g 


ii 

3t 


(44) 


then applying Fourier analysis to the record data of u, v, and p, the first 
five Fourier components of the directional spectrum can be obtained. 

e. Other Methods 

99 

Other methods, such as use of sun glitter as reported by Stilwell 
have not been fully evaluated. The 180° ambiguity is inherent. 
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APPENDIX A 

This Appendix contains a number of waveforms which resul-ted from the 
computer simulation study of the altimeter. Each series of figures is 
preceeded by a short description of the simulation parameters. 
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Figures A-l through A-4 illustrate typical simulated results of 
2 

square-law (E ) and linear (E) detector waveforms for a 50 ns rectangular 
pulse. These data correspond to the backscattered signal for the limiting 
case of a fast rise-time transmitter pulse, wide-band receiver, and 
infinite signal-to-noise ratio. 


A-2 



A- 3 














A-5 















Figures A-5 through A-9 illustrate typical simulated results of 
2 

square-law (E ) and- linear (E) detector waveforms for a 50 ns pulse of 
Gaussian shape and for noise-free reception. These data are representative 
of the waveforms for a matched- filter receiver. 
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Fig. A-5. Typical Simulated Results of Square Law and Linear 
*• Detector Waveforms for a 50 ns Gaussian Pulse, SNR = » 
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Fig. A-6. Typical Simulated Results of Square Law and Linear 
Detector Waveforms for a 50 ns Gaussian Pulse, SNR = » 




Fig. A-7, Typical Simulated Resu! 
Detector Waveforms for a 50 ni 





Fig. A-8, Typical Simulated Results 
Detector Waveforms for a 50 ns G 




Figures A-9 and A-10 are comparable to the data shown In Figure 2-3 of the 
text, except that new random numbers were "used for both signal and 
noise in Figure A-10. 
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Fig, A-10. Typical Simulated Results of Square-Law and Linear Detector Waveforms 
For a 50 ns Gaussian Pulse, SNR = 20 db (Second Computation) 






This section contains several double-delay differencing results. Figures A-ll 
and A-12 show the effect of receiver noise bandwidth and comparison of 
Figures A-ll and A-13 demonstrates the effect of new random numbers (for 
both signal and noise) on the waveforms involved in double-delay 
differencing. Additional double-delay differencer output waveforms are 
shown in Figures A-14 through A-16. Figures A-17 and A-18 show the 
double-delay differencer output for a Gaussian -shaped pulse, when the 
input is a 50 case average for SNE. =co and 30 db, respectively. 
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Fig. A-17. Double-Delay Differencer Output When 
Input is Average of 50 Waveforms, SNR = « 
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Fig. A-18. Double-Delay Differencer Output When 
Input is Average of 50 Waveforms, SNR = 20 db 
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Figures A-19 through A-26 show the results obtained by averaging 50 
individual waveforms , for various signal and noise conditions . 
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APPENDIX B 

THEORETICAL ALTITUDE MEASUREMENT ANALYSIS 

This appendix contains the results of a theoretical analysis of 
altitude errors for a double-delay processor type of altimeter. The 
transient region of the sea-return signal is modeled by the technique 
shown in Fig. B-l. 



Fig. B-l. Signal generation model. 

This signal is then combined with thermal noise, fed into a square 
law detector, and finally processed by a double-delay differencer and 
zero-crossing extractor. Exponential type autocorrelations were used in the 
analysis because of the difficulty in obtaining closed-form expressions 
with other functions such as the Gaussian. The purpose of the signal model 
(Fig. B-l) is to account for pre-detection filtering characteristics in 
the return signal. These details will be given first. 

Referring to Fig. B-l, the input autocorrelation function for the 
noise is assumed to be 


y'i'V “ s o SCt rV 0 < vb- i T • ^ 

The output of the gating function is nonstationary noise with the 
properties 


s o 6(t 1 - t 2 ) uctp U(t 2 ) 0 < t ± ,t 2 < » 


(B-2) 


B-l 






where U(*) is the unit step function. For an IF filter impulse response 
e U(t), the output correlation R yy is given by 
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The integration of (B-5) gives, 
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Further Algebraic substitution yields the result. 
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Fig. B-2. Detector model. 


Assuming that the signal is a (non-stationary) narrow-band Gaussian 
process to which the thermal noise 


R (t) = N e 

nn o 


-B |x 
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is added, the output of the square-law detector z(t) (see Fig. B-2) is 
derived as follows: 
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z(t) = [y(t) + n(t) ] 


R zz (ti,t 2 ) = E{z(t 1 ) z(t 2 )} 


2 2 
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indicates uncorrelated term. 
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. Collecting the previous results, the mean square output 'of Square Law 
Detector is 
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The signal-to-noise ratio (SNR) at detector output is. 
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Some of the above results have been given for purposes of documentation. 

The autocorrelation at the output of the double-delay differencer is readily 
derived, however, the .zero-crossing solution has not been obtained. 
Approximate results will now be discussed. 

Referring to the idealized waveforms shown in Fig. B-3, the variance 
at the instant of zero-crossing can be estimated as follows: 



Figure B-3. Differencer characteristic. 


Only two stages of the differencing operation are necessary for the discussion. 
The variance at the point zero crossing t^ is for fc 2 > t l ’ 
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where u z is the mean value, i.e.. 
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Substituting from (6) , (8) and (9) with 
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and for the matched conditions. 
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the variance at the zero crossing normalized to a unit mean value is 
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This amplitude variance can be scaled to a time variance as 
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where the mean slope is ^ rom B-3 approximately ^-/volt 

(T = pulse length, in n.s.). Using an altitude uncertainty of 

CT 

cr = -t~ a 
a 2 t 

the final altitude error estimate is 

. 0.15T 

cr = 

a /"S’ 

where n is the number of samples. 
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Appendix C 


This appendix presents the results of a brief survey' of prior work in 

estimation theory from the standpoint of the altimeter problem. The objective 

* 

of the survey is to compare optimum altitude extraction schemes 

with the ad hoc threshold, split-gate, and delay differencing techniques. 

Since the altimeter signal constitutes return from the sea-surface with 
simultaneous range and doppler spreading, the target scattering function 
o(x,f) approach is adopted. The received signal characteristics are 
considered prior to discussing optimal processing techniques. 

The basis for assuming that the scattered signal fluctuations will 
be Rayleigh distributed (at the output of a linear detector) was discussed, 
in Section III of this report. The^iysical modeling of the scattering 
process relative to a first-order convolutional model of the altimeter 
signal can be established as follows: The signal s(t) reflected from a 

small region of the sea-surface, with a two-way time delay x , will 
be a sample function of the clutter process, i.e.. 


s (t) =0 s (t— x ) (1) 

n n n 

where a , is proportional to the scattering magnitude and inversly 
proportional to range squared, within the region illuminated. The total- 
received signal during a time interval T is obtained by summing terms 
i.e. , 


s o ( t) 
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" I 

n~l 
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S(t-T ) 
n 



< 
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Since 


s(t) * <5(t - t) = 


s(t), «5(t~t -t) 


dt^=s(t-x) 


this result can be expressed as the convolution (*) of s(t) with an impulse 

train 6(t-x ), 
n 


"optimum" in the sense that some index of performance is minimized 
(or maximized) . 


C-l 



Therefore 


N 

s (t) * I<r n s(t) * <s(t~T n ) f (3] 

n=l 

which can formally be written as 

N 

s q (t) = s(t) * £ cr^ <5(t-T n ) . (4) 

c ‘ n=l 

This last form places in. evidence the impulse response nature of the sea 

return. The impulse response will be non-stationary if the range X°r time) 

dependence of theca’s is taken into account. The stationary assumption 

should certainly be valid for the altimeter problem. Equation 4 shows 

that the impulse response of the sea surface can be conceptualized as 

arising from a collection of discrete reflectors of sizes 0 and delay t . 

n J n 

This is essentially the concept used in the simulation. Use of terms such 
as impulse response of spread function to describe sea scattering are not 
precise, since these idealizations can be described in the mean only. Also, 
the impulse response cannot be rigorously defined in terms of scatter 
obtained by illuminating an infinitesimal area containing a normal incidence 
diffraction element. Equation 4 is therefore considered to consist of a 
convolution of a signal with a pseudo— impulse response. 

For the altimeter problem, the clutter impulse response from (4) is 

N 

h c (t) = l a n 4(t“T n ). (5) 

n=l 

Ocean surface correlation has been neglected thus far in the discussion. 
Arguments given in Section IV indicate that the sea return correlation time 
at the altimeter can be less than one nanosecond. Therefore, for the 
present purposes the delta function concept will be valid for time-increments 
of one nanosecond or greater, or signal bandwidths less than about 1000 MHz, 
for normal-incidence altimeter geometry. On the above basis, samples of the 
sea-surface impulse response can be generated by selecting random numbers. 
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Assuming that the clutter is "frozen" during the period of the r-f pulse 
(rectangular of length T,) the scattered signal would be of the form 


h (t - T) * [U(t) - U(t - X)] 
c 


where U(t) is a step function. 

The pulse response for a particular sample function of h^ thus consists 
of a transient period (the integral of the impulse response) followed by a 
steady state period in which the output is a moving average over the time 
interval T. The sample functions given in Fig. 2-1 are examples of 
detected sea-return for rectangular pulse illumination. 

This model of sea-return also leads to the Rayleigh description of 
radar return. The scattered signal will be a sum ofphasors. 


a 

n 



This expression, which consists of a random collection of phasors, leads 
to the well known Rayleigh probability distribution. Since the step 
response is the integral of the impulse response, the Rayleigh model 
can be shown to lead to a time-varying Gaussian envelope distribution of 

the for " 1 I 2 1 

p(v> “■> |S?t| ‘i 1 

'/,’h fire Y is the process variance. P(v) will be - a Wiener Levy process^ 
in which the variance increases linearly with time. For the above 
model the spectrum (the Fourier transform of h c (t)) will be nearly 
constant up to some frequency on tbe order of the reciprocal of the correlation 
time (<1 ns). The bandwidth requirements for a signal which accurately 
probes the sea-return impulse response can thus be estimated. Visualizing 
the above described spectrum as one corresponding to a linear system, it 
is obvious that the exploring signal must contain frequency components 
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that extend throughout the transform of h c (t) • Therefore a signal that 
is designed to make complete measurements of the impulse response 
sample functions must contain frequency components in the gigahertz 
range. The measurement of this quantity is of concern in the radar 
altimeter design. As discussed in Section II, one of the problems in 
design of a system which attempts to estimate the mean value of a fluctuating 
signal, is the problem of obtaining an adequate number of independent 
samples. For a transmitted pulse length in the order of 50 nanoseconds, 
errors arising from self noise appear to be comparable to thermal noise errors. 
Therefore design of a precision altimeter may greatly benefit from 
improved knowledge of the characteristics of the received signal. Conventional 
(non-adaptive) matched filter IF bandwidth criteria may not be appropriate 
since, for unequal errors in the two sources, an optimum bandwidth would be 
expected to exist. 

With these qualitative remarks on signal processing considerations, 
the subjects of detection and parameter estimation will be considered. The 
parameter estimation problem is simplified by considering signals 
existing in the IF, since the statistics are time-varying Gaussian at 
this point. Otherwise, non-Gaussian statistics must be considered. 

Neglecting, for the time being, the non -deter min istic nature of the 
scattered signal, the optimum processor for measuring range (time of arrival) 
in the presence of additive Gaussian noise will be considered. For a 
partially coherent system the processor shown in Fig. C-l is optimum. This 
processor consists of quadrature detectors, matched filters, and a bank 
of correlators ^ ^ . 

The limitation of the above theory is that it does not take into account 
the random nature of the received signal. The fluctuating altimeter 
signal is strongly akin to multipath signals encountered in radar 
astronomy and scatter communications. The solution to this problem, due to 
Kailath^ , is diagrammed in Fig. G-2. The estimation filter H shown in 
Fig. C-2 is derived from the equation 
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H = (1+A A ) 
n sy 


where A is the noise covariance and A is the covariance of the channel 
n sy 

output y(t) . 

The scattering function 0 " (x) approach developed by Price and 

Green and discussed by Evans and Hagfors 7 will be reviewed next. This 

approach is based on the likelihood ratio of the form 

.T 

Q 

|| W o (t) W o (t’) <j> g (t,t') dt dt' 

-T 

o 

2 


Where T 

o 

W 


o 


4 > 


g 


is the observation time, 

is a sample observation of the received signal, and 
is the echo correlation function. 


This last quantity is related to the scattering function through the 
relationship 


<j> g (t,t’) = ~ Re [exp tjW Q (t - t')] j X(t - x) X*(t’- x) L y (t - t’, x) dx] 

—CO 

in which X(‘) is the complex transmitted envelope and L is given by 


Iy(At,x) 


a (w,x) e 


jwAt 


dw 


The solution to this set of equations provides what is termed the 
"practically optimum" processor . The optimum range estimation technique, 
for a target spread in both delay and Doppler, consists of a bank of 
processors spaced in range and the receiver chooses the processor exhibiting 
the largest output. This filter shown in Fig. C-3 is termed a weighted 
radiometer, in contrast to a matched filter envelope detector. The first 
filter shown in this figure is one matched to the transmission replica 
s(t). This is followed by a squarer (envelope detector) and a filter matched 
to the mean target delay characteristics . The impulse response of the latter 
is 


07 



h(t) = H 2 (x^-t) 

2 

where H (t) is the square of the filter* matched to the scattering function 
cr (u o> t). Thus for fluctuations small relative to the spectrum of s(t), 
the impulse response is essentially o(t), This form of processor is shown 
in Ref. 7 to be opt imum if 1) a Gaussian pulse is used, 2) a(T 3 w) is 
factorable in rand m, and 3) if the observations extend over all time. 
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Fig. C-3. 


The implementation of such a processor in the altimeter does not 
appear promising. The only novel feature is the video filter matched to the 
target time-profile. Based on our previous discussion, the filter would 
have an impulse response of the type shown in Fig. 4-C. 



In practice, realizability would require some time delay. However, its 
general form appears to be an "integrate and dump" circuit, which is 
similar to the split-gate concept studied on the altimeter program. 
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This discussion indicates that the synthesis of the above type of 
optimum processors will not yield substantial improvements over the present 
ad hoc techniques. The areas of waveform synthesis and adaptive filtering 
should be investigated. 
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APPENDIX D 


Analysis of Ocean Surface Effects on. the Received Waveform 

In order to account for the rough sea surface, it is necessary to 
compute the density of stationary points per unit area as a function 
of height above the mean sea surface. Rice‘S has examined the distribu- 
tion of local maxima and minima for a bandlimited signal, which corres- 

2 

ponds to a one-dimensional sea. Longuet-Higgins examined the density 
of specular points with height above the sea surface. Making use of 
the latter work and the fact that the density of stationary points 
(specular points for normal incidence illumination) is four times the 
density of the maximum points, the distribution of stationary points per 
unit area may be established as a function of height above mean sea level 
(msl) . 

Formally, the required probability is given by 


4dz 


*OrC 

—co —oi 


z z expl^T^T [M^z 2 + 2M_. z + z ML _ 
XX yy r (2 M| 11 14 xx yy 15 


+ 2z z M /c + z 2 M. . + z 2 dz dz 

xx yy 45 xx 44 yy 55 j xx yy 


where dp is the probability that a maximum exists between z and z + dz. 

z is the height of the surface above msl (a function of x and y) , the 

position coordinates on the mean surface, and 

*2 

_ a z 

Z XX 2 

ax 


z yy a 2 

dz 

and M, . are elements of the correlation matrix for the random variables 
ij 2 2 

z, ~, 7 ^, — f - , |m| is the determinant of the correlation matrix. 

The elements of the correlation matrix can be related to derivatives 
of the correlation function of the surface (and hence to the energy 
spectrum) . For example: 
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This approach cannot be pursued because the moments involving k and 

4 

k are not known to sufficient accuracy. For presently available oceano- 
graphic spectra, the last three integrals are infinite. 
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